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Abstract. We diagonalize the anti- ferroelectric XXZ-Hamiltonian directly 
in the thermodynamic limit, where the model becomes invariant under the 

action of t/ q (sl(2)). Our method is based on the representation theory of quan- 
tum affine algebras, the related vertex operators and KZ equation, and thereby 
bypasses the usual process of starting from a finite lattice, taking the thermo- 
dynamic limit and filling the Dirac sea. From recent results on the algebraic 
structure of the corner transfer matrix of the model, we obtain the vacuum vec- 
tor of the Hamiltonian. The rest of the eigenvectors are obtained by applying 
the vertex operators, which act as particle creation operators in the space of 
eigenvectors. 

We check the agreement of our results with those obtained using the Bethe 
Ansatz in a number of cases, and with others obtained in the scaling limit — 
the sw(2)-invariant Thirring model. 
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§0. Introduction 

0.1 A diagonalization scheme In this paper we give a new scheme for diagonalizing 
the 1-dimensional XXZ spin chain 

1 oo 

Hxxz = — 2 E K+W + + A ^+i^)> (O- 1 ) 

/c— — oo 

for A < —1, directly in the thermodynamic limit, using the representation theory of the 
quantum affine algebra U q (sl(2)): we consider the infinite tensor product 

W = ---®C 2 ®C 2 ®C 2 ®C 2 (0.2) 

on which Hxxz formally acts, rather than starting with a finite product and subsequently 
taking the thermodynamic limit. The reason is that the model is [/ 9 (s[(2))-symmetric 
only in the thermodynamic limit, and the presence of that symmetry is central to our 
approach. On the other hand, working directly in that limit makes it difficult to give 
rigorous analytic proofs, and a number of our statements will be given as conjectures, 
supported by explicit calculations. 

Our scheme is a direct descendant of the recent works [1-4], and Smirnov's picture of 
the stt(2)-invariant Thirring model [5], which is the continuum limit of the XXZ model. 
(See also [6,7].) More broadly, it originates in many of the developments that took place 
in the past two decades in the field of integrable models in quantum field theory and 
statistical mechanics. It is remarkable that in the context of this simple off-critical model, 
one can recognize many of these developments, and for the rest of this introduction, we 
wish to recall certain aspects that are relevant to our work. 

The peculiarities of this branch of science are twofold: mathematically, its signifi- 
cance is that the systems we are dealing with have infinite degrees of freedom (unlike 
those in conventional mathematics, in which we manipulate finite degrees of freedom), 
while physically, it is the fact that they are integrable (unlike most of the physical prob- 
lems) because of their infinite, and sometimes hidden symmetries. The symmetries serve 
as the magical word in Arabian Nights that opens the door between the two worlds of in- 
finite and finite degrees of freedom. Thus, through the study of the integrable models we 
have been discovering many unexpected links between different branches of mathematics 
such as representation theory, differential equations, combinatorics, topology, algebraic 
geometry, and so on. 

Let us discuss specifically the integrability and the symmetry of lattice systems. For 
two-dimensional lattice systems, the integrability is commonly understood as the exis- 
tence of a family of commuting transfer matrices. In the language of one-dimensional 
quantum chains, it is understood as the existence of infinitely many mutually commuting 
conservation laws. Let us call them the infinite abelian symmetries. Because they are 
commuting we can simultaneously diagonalize them, and in most cases the spaces of 
common eigenvectors are finite-dimensional. In other words, the infinite abelian symme- 
tries reduce the degeneracies of the spectrum from infinite to finite. This is the reason 
for the solvability. But it is not the whole story. In this introduction, wc will discuss 
other types of symmetries that we will refer to as the non-abelian symmetries and the 
dynamical symmetries. 

0.2 The Bethe Ansatz The Bethe Ansatz was invented in order to solve the XXX 
model, then it was applied to the XXZ and many other models. Let us discuss the XXX 
model specifically. The Bethe Ansatz reduces the diagonalization of the Hamiltonian, a 
2 N x 2 N (huge!) matrix where N is the size of the system, to a system of m coupled 
algebraic equations, where m is less than or equal to N/2. For a finite but large N, this 
system is very complicated, and far from "solvable" . (It is not completely settled even for 
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m = 2 [8].) On the other hand, and surprisingly, in the thermodynamic limit (N — ► oo) 
the system of algebraic equations changes into a system of linear integral equations that 
can be solved. In this way, the ground state and the low-lying excitations have been 
extensively studied for the XXX model, and by a similar manner for many other models. 

In [9,10] the anti- ferroelectric XXX and XXZ chains are studied in the thermody- 
namic limit N = oo and at the physical temperature T = 0. We also deal with the 
same problem in this paper, but by a completely different method. Our diagonalization 
scheme is totally independent of the Bethe Ansatz, and goes further: it makes it possible 
to describe all the eigenvectors by using the powerful tools of representation theory. 

0.3 Non-abelian symmeries and particle picture In a remarkable paper [11], Faddeev 
and Takhtajan conjectured the structure of the eigenvectors of the XXX Hamiltonian. 
Using the Bethe Ansatz, they conclude that all the excitations are composed of elemen- 
tary ones, called particles, and that the physical space is 



(0.3) 

symm 

where [ ] sy mm stands for an appropriate symmetrization. In [11] the summation is only 
for even n. In our work, we consider both even and odd n. 

The issue in [11] is as follows: The XXX model possesses an sl(2) symmetry. The 
Bethe Ansatz produces only the highest weight vectors of this st(2). Therefore, spin-1/2 
two-particle states (i.e., the case n = 2, C 2 <8>C 2 ), which were correctly observed in [11] in 
the thermodynamic limit, appear in the finite chain to decompose into spin-1 plus spin-0 
particles (i.e., C 3 © C). In fact, the energy levels of the triplet C 3 and the singlet C are 
different for finite lattice and become equal only in the thermodynamic limit. This gives 
an example of how a system can achieve extra and big symmetries in the infinite volume 
limit, and shows one of the advantages of a direct diagonalization of the Hamiltonian in 
the infinite lattice. 

Let us briefly consider the particle picture in the XXZ case. By particle picture 
we mean the structure of common eigenspaces as irreducible modules of the algebra 
of non-abelian symmetries. Here, we mean by non-abelian symmetries operators that 
are commuting with the Hamiltonian. As opposed to the abelian symmetries, they 
are non-abelian, so they change eigenvectors of the Hamiltonian without changing their 
eigenvalues. The finite XXZ chain has no s[(2) symmetry, but as Pasquier and Saleur 
[12] pointed out, if we add certain boundary terms, then the XXZ Hamiltonian acquires 
the C/g(sI(2)) symmetry. But the representations are still highly reducible. 

In our approach, instead of adding boundary terms to the finite chain, we consider 

the infinite chain from the outset, so that the J7 g (sl(2)) symmetry (which is much larger 

than [7 9 (s[(2))) is manifest. As J7 g (s[(2))-module the common eigenspaces of the XXZ 

and the higher Hamiltonians become the irreducible tensor products £g)™C 2 parametrized 
by n quasi-momentum variables. This is exactly the point of Faddeev and Takhtajan (in 
the XXX-c&se), though they did not explain it in these words. 

0.4 Dynamical symmetries and creation and annihilation operators Since Onsager's 
solution of the Ising model, a number of alternative methods of solving that model were 
developed. A method using the infinite abelian symmetries is given, e.g., in Baxter's 
book [13]. A major difference, when compared with the anti-ferroelectric XXZ model, is 
that the Ising model has non-degenerate common eigenvectors of the commuting transfer 
matrices. In other words, the Ising model has no non-abelian symmetries which com- 
mute with the Hamiltonian. However we may also consider symmetries in a broader 
sense: those which do not commute with the Hamiltonian but map eigenvectors to other 
eigenvectors with different eigenvalues (let us call them the dynamical symmetries) . The 
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Ising model has such symmetries: they are generated by the creation and annihilation 
operators of free fermions. 

A natural question arises: what are the dynamical symmetries of the XXZ model? 
Our answer is that the vertex operators for C/ g (s[(2)) give an appropriate mathematical 
language for them. We will explain this statement later. 

Of course, the lesson of the Ising model is not restricted to the fact that the space 
of particles is the Fock space T of free fermions. The most striking fact observed in 
the Ising model is that the spin operator s n is completely characterized by its adjoint 
action on T, i.e., Ads„ G End (J-), or on even smaller linear subspacc spanned by the 
free fermions. (2 N x 2^ reduces to 2N x 2N\) This was the key (or the magical word) to 
the world of the monodromy preserving deformations, the Painleve transcendents, etc. 
[14,15]. 

Around the same time as the above developments in the Ising model, the theory 
of the S'-matrix was developed by Zamolodchikov and Zamolodchikov, on the basis of 
the fact that if a field theory has infinite abelian symmetries, then it has a factorized 
S-matrix. They defined the algebra of the creation and the annihilation operators as 
the key to the bootstrap program for the determination of the S-matriccs. The creation 
and annihilation operators, which we mentioned above as generators of the dynamical 
symmetries of the XXZ Hamiltonian, give a lattice realization of the Zamolodchikov 
algebra [1]. 

0.5 QISM The quantum inverse scattering method (QISM), initiated in Leningrad 
when the city was called by that name, is a large-scale project to understand intcgrable 
quantum systems as a whole. Among its many achievements, we recall two that are 
related to this work. One is, of course, the discovery of quantum groups, or the q- 
analogue of the universal enveloping algebras, by Drinfeld and Jimbo. The other is 

Smirnov's work, to which we come later. The reason for the C/ g (sI(2)) symmetry of the 
XXZ Hamiltonian is simply that the Boltzmann weight of the six-vertex model is the 
i?-matrix for the two-dimensional representation of U q (sl(2)) (i.e., C 2 in (0.1)). Also, 
the existence of the infinite dimensional abelian symmetries is a corollary of this fact. 
The discovery of quantum groups came as a harvest of QISM, not conversely. All the 
developments in the earlier days of QISM were made in the absence of quantum groups. 
Therefore, they are lacking in understanding of the true nature of the symmetries. The 
QISM project should be further pushed forward on the basis of the developments in the 
theory of quantum groups. 

In the early days of QISM, Smirnov invented a bootstrap program for computing 
the matrix elements of local fields starting from a factorized S'-matrix [5,16]. He found 
a system of difference equations which ensures the locality of the field operators, and 
obtained the form factors (i.e., the matrix elements of the off-shell currents). This 
is a deep result. We are only beginning to understand its true meaning. Frenkel and 
Reshctikhin obtained the q-deformed Knizhnik-Zamolodchikov equation, and established 
that the n-point correlation functions of the vertex operators satisfy this system of q- 
difference equations [1]. Smirnov's equation was the double Yangian version of that. 

As we already mentioned, the creation and the annihilation operators of the XXZ 
model are given in terms of the vertex operators of J7<j(s[(2)). (They are not equal. 
See section 7 for a detailed discussion on this point.) The latter are exactly the same 
as special cases of those considered in [1]. So, their n-point functions satisfy the g-KZ 
equation. This is the key to the calculation of the energy and the momentum, and the 
commutation relations of the creation and annihilation operators. 

0. 6 Conformal field theory We have considered the symmetries of the Ising model and 
the XXZ model. The importance of such a viewpoint was established in the monumental 
work of Belavin, Polyakov and Zamolodchikov [17]. To formulate the symmetry picture 
of conformal field theory in a way that is suitable to our purposes, we will consider the 
st(2) Wess-Zumino model, following Tsuchiya and Kanie [18]: 
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(i) The space of states of this model is a direct sum H of the irreducible highest weight 

representations of st(2). 

(ii) The Hamiltonian of the model is the Virasoro generator Lo in the Sugawara form. 

It allows the sl(2) symmetry. The total sl(2) acts on H as dynamical symmetries. 

(iii) The local fields of this model are given by the vertex operators: 

$(z) :V Z ®H->H. (0.4) 

Here V z is a level representation of sl(2) depending on a parameter z, and $ 
intertwines the representations on both sides. 

(iv) The vacuum-expectation values of the vertex operators are characterized by the 
Knizhnik-Zamolodchikov equations, which is a system of linear holonomic differen- 
tial equations. 

Now we can give the definition of the vertex operators for C/ 9 (s[(2)). They are 
exactly (0.4) in the context of the quantum affine Lie algebra £/ g (s[(2)). Note, however, 
that there is a difference between the interpretations of the spectral parameter z in (0.4), 
in the WZ-model and in our XXZ model. In the former, it is the coordinate variable, 
and in the latter, it is the momentum variable. The mathematical extension from q = 1 
to q ^ 1 of the vertex operators has a different physical content. (This is already pointed 
out by Smirnov in a different context [5].) 

Conformal field theory has had many achievements in the subject of integrable mod- 
els. In fact, we should say that it has changed the definition of the game we are playing. 
However, it has not made progress in all of the subjects studied in earlier days. The 
reason is obvious: its basic principle, the conformal symmetry, is valid only in massless 
theories. Zamolodchikov made the first attempt to attack the massive theories from 
the conformal viewpoint. He proposed to define integrable deformations of conformal 
field theories by the existence of infinite abelian symmetries embedded in the conformal 
symmetry [19,20]. This motivated a reconsideration of the factorized S'-matrix theory 
[21]. 

0. 7 CTM and beyond Let us return to the XXZ model. The point we wanted to 
make is the following: The deformation parameter q appears in the XXZ model as the 
anisotropy parameter A = (q + q~ 1 )/2. This is the departure from criticality, or simply 
the mass. So, there is more than a good reason to expect that the quantum affine algebras 
are relevant in massive integrable models. In fact, many of the recent developments are 
related to this point. For the continuum theory, we refer the reader to Smirnov's paper 
[5]. 

Now let us come to one of the masterpieces of Baxter, the CTM (corner transfer 
matrix). In [22], after saying that there is no Ising-like reduction from 2^ x 2 N to 
2N x 2N in the six or eight vertex models, he wrote 

A rather ambitious hope is that by examining the CTM's we may stumble on 
such a group, that the solution of the models may thereby be simplified . . . 

A few years later, he succeeded in computing the 1-point functions of the hard-hexagon 
model and its generalizations, without giving an answer to the original question of his 
earlier paper. Nevertheless he has given us the magical word — "Open Sesame", by 
inventing the CTM method. When the cave-door opened, there appeared the characters 
of the affine Lie algebras [23,24], the theory of crystals [25,26], the irreducible highest 
weight representations [3], the q- vertex operators [1,4] — the passage to the represen- 
tation theory of the affine quantum groups. And in every case, the use of an infinite 
system, already assumed in Baxter's original work, is an essential ingredient. 

In this paper, we proceed further along this route. We propose that the mathemat- 
ical picture that bridges (0.2) and (0.3) is 



End c (V(A )) = V{K a ) ® V (A )* 



(0.5), 
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(in the even particles sector; the odd particles sector is Horn c (V(A ), V(Ai)) = V(Ai) ® 
V(Ao)*), where the V(Ai) are the level 1 highest weight representations of i7 g (s[(2)). 
By this, the world of the infinite degrees of freedom (i.e., the infinite tensor product) 
opens to the world of the finite degrees of freedom (i.e., the representation theory of 

?7q(s[(2))). The vacuum, the lowest eigenvector of the Hamiltonian, is the identity 
operator in End c (V(Ao)) , or the canonical element in V(Ao) (S> V(Ao)*. The particle, 
as given in (0.3), are created by the vertex operators acting on the left half of V(Ai) ® 
V(A )* (i = 0,1). Thus, the space T (0.3) of the eigenvectors of the Hamiltonian 
(0.1), that are created by the creation operators upon the vacuum, lies in the level zero 

U q (sl(2)) -module (0.5). 

Let us summarize our discussion on the symmetries of integrable models in this 
introduction. 



Symmetries Ising XXZ Thirring WZ 

Infinite None 
Conservation Laws 

None U' q (s[{2)) The s[(2)-Yangian su{2) 



Abelian Commuting 
Symmetries Transfer Matrices 



Non-abelian 
Symmetries 



Dynamical 
Symmetries 

Space of States 



Free 
Fcrmions 

T 



Vertex 
Operators 

T C End(W) = H®H* 



Local Fields Clifford group 



Virasoro algebra 
su{2) 

U 

Vertex Operators 



0.8 Plan of the paper The text is organized as follows. In Sect. 1 we introduce the 
quantum affine algebra U q (g) as non-abelian symmetries of the XXZ spin chain. In Sect. 
2 we describe the embedding of the highest weight modules V(Ai) into the half-infinite 
tensor product space • • • V ® V ® V or V ® V <8> V <8> ■ ■ ■ given by iterating the vertex 
operators. Examining the perturbation expansion at q = we observe that, choosing 
the scalar multiple of vertex operators correctly the scries defining the embedding would 
become finite, term by term in powers of q. To ascertain the existence of such a nor- 
malization factor we need to know the aysmptotics of the n point functions of vertex 
operators (cf. Sect. 6.8 below) as n — > oo. Though the problem is of its own interest, it 
is beyond the scope of the present paper. In Sect. 3 we study the decomposition of the 
level module V(A ) <g> V(A )* at q = 0. We show that its crystal naturally decomposes 
to 'n-particle sectors', each of which can be identified as a certain anti-symmetrization 
of the affine crystal Aff(B)®". We study the picture for nonzero q in the next sections. 
In Sect. 4 we introduce the vacuum vectors as canonical elements of V(Ai) ® V(Aj)*. 
In Sect. 5 we formulate the 'particle picture' by utilizing a similar but different kind of 
vertex operators from those used in the embeddings. We argue that in this setting the 
computation of the one-point function (in the sense of [23] , not the staggered polarization 
[27]) can be done trivially. 

Sections 1-5 are devoted to a presentation of our ideas on the problem. In the 
next two sections we reformulate the problem in the language of representation theory. 
After reviewing the g-deformed vertex operators following [1] and [4] we give formulas 
for their two point functions and their commutation relations. The case of general n 
point functions is discussed in Sect. 6.8. In Sect. 7 we study the vacuum, creation and 
annihilation operators which arc defined in purely representation-theoretic terms. We 
prove that the vacuum vector has the correct invariance with respect to the translation 
and energy operators. We then derive the energy-momentum of the creation-annihilation 
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operators using the formulas for two point functions, and show that they coincide with 
the known results obtained by the Bcthc Ansatz method. The Sect. 8 is devoted to a 
summary and open problems. 

In the appendices we collect some data concerning the global crystal base and vertex 
operators. We give tables for the first few terms of the actions with respect to the 
global base of U q (g) on V(Aq) (Appendix 1), of the vertex operators (Appendix 2), the 
embedding into the infinite tensor product (Appendix 3) and the images of the vacuum, 
one- and two-particles (Appendix 4). In Appendix 5 we study the q — > limit of the 
Bcthe vectors and compare the results with the vertex operator computations. 

§1. Quantum Affine Symmetries of the XXZ model 

Let V ~ C 2 be a 2-dimensional vector space, and let <r x , a v , a z be the Pauli matrices 
acting on V. We consider the infinite tensor product W (0.2) and the XXZ Hamiltonian 
Hxxz that formally acts on W . The action of Hxxz is formal in the sense that it is a 
priori divergent, since we are working directly in the thermodynamic limit, and requires 
renormalization. We number the tensor components by k € Z. The limit k — > oo 
is to the left, and k — > — oo to the right. Our aim is to diagonalizc the Hamiltonian 
Hxxz, using the representation theory of quantum affine algebras. The Bethe Ansatz 
provides us with a method for diagonalizing such Hamiltonians. In this method we start 
from a finite (periodic or non-periodic) chain, find the eigenvalues and the eigenvectors 
of the Hamiltonian in a certain form, and take the thermodynamic limit at the end. 
In this paper we propose another method that diagonalizcs Hxxz directly on W . The 
applicability of our method, at the moment, is limited to the anti-ferroelectric regime A < 
— 1. The Bcthc Ansatz method in this regime meets with the difficulty of starting from 
the wrong vacuum, and filling the Dirac sea. Our method is free from this complication. 

The main ingredient of our method is the quantum affine symmetries of U q (sl(2)). 

Consider the action of U' q (s{{2)) on V; 



tt : ^(st(2)) 



End(V) 



given by 




We set 




We relate q to A by 



A = i^, -K ff <0. 
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By using the infinite co-multiplication A(°°), we define the action of U q {si{2)) formally 
on W: 

A (oo) ( ei ) = ^•••®t i ®e i ®l---, (1.1) 
kez 

A(°°)(/ i ) = ^---®i®/ i ®tr 1 ---, (i.2) 

kez 

A( tx> )(t i ) = ---(g)t i ®t i (8)t i "-. (1.3) 
One can check the commutativity 

[H xxz ,U' q (si(2))}=0 

as discussed in [3]. The major advantage of working in the infinite lattice limit directly is 
having this infinite symmetry, which inevitably breaks down on a finite lattice. On a finite 
lattice, one can modify Hxxz in such a way that it possesses the C/ 9 (s[(2))-symmetries 

[12], but not the U q (sl(2)) -symmetries. (Otherwise, the degeneracy of the spectrum on 
the finite chain would be different.) The difference of the size of the symmetry algebra is 
crucial. For instance, we can identify the vacuum vector in W as the unique U' q {s[{2i))- 
singlet (a vector which generates a 1-dimensional t^(s[(2))-module.) On the other hand, 
singlets for U q (sl(2)) are many. 

The full algebra U q (st(2)) is obtained by adding q d to U' q {s\{2)). Following [3] let 
us identify d with (Hctm — S)/2 where 

hctm = --^ J2 K° x k+i° x k + °l + l<% + ^t+lO (1-4) 

1 kez 

and 2S — J2 a k ^ s tne total s P m operator. This identification is justified by explicit 
computations and by checking that 

[d,e l ] = S^, [d, fc] = -6 M fi, [d,t»] = 0. 

Consider the shift T of W which induces the outer automorphism of the operator 
algebra generated by cr£, a\, cr|; 

T is a shift to the right by one lattice unit. Because of (1.4), the Hamiltonian (0.1) can 
be written as 

—-^Hxxz = T ■ d - T^ 1 — d. (1.5) 
1 — q z 

This is the key to our diagonalization procedure: the point is that the derivation d 
[3] , as well as T can be expressed in the language of representation theory, thus the XXZ 
Hamiltonian can be expressed in terms of mathematically well-defined objects. 

In fact, (1.5) allows us to bypass explicit references to the higher Hamiltonians, 
obtained by taking higher derivatives of the row-to-row transfer matrix of the six- vertex 
model. The point is that all higher Hamiltonians can be generated from i.e., Hxxz, by 
taking commutators with d recursively [28]. This means that the essential feature of 
the integrability (i.e., the existence of the infinite conservation laws), in this model, is 
incorporated in d. 
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Let V(Ai) (i = 0,1) be the level 1 irreducible highest weight [/ g (s[(2))-module with 
highest weight Aj. We give a conjecture on an embedding of V(Ai) into the half infinite 
tensor product 

Wi = ■ ■ • ® V ® V 

by using the vertex operators 

$ : V(A t ) — > V(Ai_i) ® V, *(«)=*+(«)<»«+ + $_(«)<»«_, (2.1) 
that satisfies 

$(ot) = A(a;)$(w) for all a; G U q (sl(2)) and w e V(Aj). 

Precisely speaking, we need a completion of V"(Ai_i) <g> V, which we will elaborate in 

Sect. 6. Also, we do not define Wi (but see the definition w[ 1 ^ below). The subscript I 
in Wi stands for left. It means the left-half- infinite tensor product. 
The idea is to consider the composition of the vertex operators; 

V(A )^V(A 1 ) ® V^V(Ao) ®V®V 

*?^V(Kx)®V®V®V — » >W l (2.2) 

Our conjecture is that with a proper normalization of (2.1), the composition (2.2) con- 
verges to a map 



i : V{Ai) — » W; 

satisfying 

t(ara) = A(°°)(a;) t (w), 

and that if v is a weight vector of V(AA, then t(v) is an eigenvector of Hctm defined 
by (1.4) with a certain renormalization (see [3]). 

Let vf^ be the set of paths which parametrise the affine crystal B{Ai). (See 

[2,29,30]) We use the convention that the colors i = 0, 1 arc modulo 2 (e.g., — 

v[ l+2) ). A pathp= (p(k)) k>1 = (•■•p(3)p(2)p(l)) G "P, W satisfies p(k) = (±) (often 

identified with ±1), and p(2k~+ i) = (+), p(2k + i + 1) = (-) for fe > 0. We denote by 
\p) the vector in Wi corresponding to p; 

\p) = ■■■®v p {3) ®Vp(2) ®v p(1) . 

The following are called the ground-state-paths; 

Po = {■■■+ - + -), Pi = ( + - +)• 

We consider the set of formal infinite linear combinations of paths with coefficients 
in Q(q), 

W® = { E c(p)|p>;c(p)eQ(«)} 
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Let {G(p);p £ V^} be the upper global base of Kashiwara [25] for V(Ai). Define 
c ± (r,p) (r.pePf'upW) by 

*±(G(p)) = ^c±(r,p)G{r) ®v±. 

r 

If p € P ; (<) , then c±(r,p) = unless r £ and s(p) = s(r) ± \. Here s : 7> (0) U 

— > \Tj is the spin of paths defined by s(p) = \ lim^oo Y^ 2 j=i P(j)- 
The coefficients c± (r, p) of the vertex operator satisfy the following (see [4] ) : 

c_(pi,p ) = c+(p ,Pi) = 1, (2.3) 
c±(r,p) = 1 mod qA if p{\) = ±, p(fc + 1) = r(fc) Vfc (2.4) 
= mod otherwise, (2-5) 

)J{r e P; (0) U | c±(r,p) ^ mod q^A} < oo 

for all p e -P ; (0) U V{1\ and ffeN. (2.6) 

Here A = {f e Q(q); q has no pole at q = 0}. (2.3) is the normalization and (2.4), (2.5), 
(2.6) are the key properties which mean the compatibility between the vertex operator 
and the crystal base. 
Define 

uW(a,p) = (G*(a n+1 )\$ a{n) o • • • o $ o(1) |G(p)>. (2.7) 

Here a = (a n +i, a(n), . . . , a(l)) and p are paths in "P^\ and {G*(p)} is the dual base to 
{G(p)}. Because of (2.4), (2.5), (2.6), (2.7) is convergent in Q[[q}]. 

Conjecture. 

(i) There exists a limit 

Uj( U H a , P) r,rr n 

w(a,p) = lim V eZ g . 

(iij w(o,p) £ Q(g). 
(iii) Setting 

< G (P)) = E w ( a 'P)l«)' ( 2 - 8 ) 



we have 



(iv) 



<G(p))\ q ^ = \p). (2.9) 



i(xG(p)) = ^ W (a,p)A(° c )(x)| a ) (2.10) 
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for x G U' q (sl{2)) and p G -P, (i) . 

The statements of the conjecture are independent of the normalization of the vertex 
operator (2.1). With the normalization (2.3) we have (see 6.8 and Appendix 3) 

Mm ^(")fe+„,K) 1/n = 1 + 9 4 - <f + 9 8 mod q w . 

n — >oo 

Therefore, we also conjecture that this is actually convergent in Z[[g]]. We have no 
reasonable conjecture what the limit is. 

In [3] an embedding V(A ) V ® V ® ■ ■ ■ is constructed by using a different 
coproduct; 

A_(e 4 ) = e i <g>tr 1 + l<8>ei, 
A_(/j) = fi®l + ti®fu 
A_(ti) = *i (8) *». 

To compare these two different embeddings, let us flip right and left in the notation of 
[3]. (Namely we change the embedding of [3] into the form V(Aq) — ► ■ • • <g) V <S> V.) If 

we denote the infinite coproduct in [3] after the flip by A^?^, we have 

A^(/ i ) = gA(~)(t i / i ), 

Therefore these two embeddings differ only by certain power of q. In Appendices 3,4 we 
compute c±(r,p) and u>(a,p) for several cases. 

Let us discuss on the similarity and the difference of [3] and the present paper. [3] 
studied the CTM Hamiltonian (1.4) on Wi and identified its eigenvectors with the weight 
vectors of V(A ) embedded in W\. In this paper we study the XXZ Hamiltonian (0.1) 

on W and identify its eigenvectors with the weight vectors in certain [/^(s[(2)) -modules 
embedded in W. In both cases the main aim is to give an equivalent mathematical 
picture to the physical content of the problem, which is given by the CTM Hamiltonian 
or the XXZ Hamiltonian. The mathematical picture presented in [3] is V(Ao), the 

level 1 irreducible highest weight U q (sl(2)) -module. The method employed is the q- 
perturbation. An obvious implication of [3] to the present problem is to consider 

V(A Q )®V(A Q )* a (2.11) 

as a mathematical model for the infinite product W. A similar analysis applied to 

W r = V ® V <g> • • • 

gives the level —1 irreducible lowest weight l/ g (s[(2))-module. The right half of (2.11) is 
the most appropriate realization of that. Here V(Aq)* = Homq( 9 ) (V(A ),Q(q)) endowed 
with the natural right U q (sl(2)) -action p^; if / G V(A n )* then (pn(x)f)(v) — f(xv) for 
x e U q (si(2)). Then V"(A )* a is the left U q (sl(2)) -module with the left U q (sl(2)) action 
PL, a in terms of the antipode a (see [2], Sect. 5); 

(pL,a( x )f)( v ) = f( a ( x ) v )- 



12 



B. Davies et al. 



We will discuss why we choose the antipode to construct the left action a bit later. Now 
we want to discuss a big difference between [3] and the present paper. The basic tool in 
[3] was the g-series expansion. Underlying this were two things; 

(i) The success of the theory of the afHne crystals. ([2,29,30]) 

(ii) The CTM magic, i.e., the discreteness of the spectrum of the CTM Hamiltonian. 
([13]) 

In this paper, we further exploit (i) to analyse the q — structure. Note that the 
XXZ Hamiltonian is already diagonal at q — 0, i.e., A = — oo. The eigenvectors are 
nothing more than paths p = (p(k))kez with appropriate boundary conditions. The 
notion of crystal gives a nice combinatorial structure to the set of paths. This is de- 
scribed in detail in section 3. The discreteness of the spectrum is no longer true for the 
XXZ Hamiltonian. It means the breakdown of the g-expansion as a tool for finding the 
eigenvectors. The continuum spectrum at q ^ degenerates to the discrete spectrum 
at q = 0. Therefore we do not know a priori the correct form of the eigenvectors at 
q = with which we should start our expansions. In [31] the expansion of the vacuum 
vector (i.e., the lowest eigenvector of Hxxz) was discussed. This was possible because 
the lowest eigenvalue is not degenerate even at q = ( or more precisely, it is doubly 
degenerate, but the two eigenvectors generate two orthogonal sectors). So we could start 
with 

\vac)\ g=Q = (... + - + -...). 

Similarly it is not difficult to tell the correct limiting form of the 1 particle state because 
the 1 particle states must be eigenvectors of T 2 (the shift of two lattice units). Therefore 
we can start with 

|«.+>U = E("-- ++ - 

k 

i^-)l g =o = E(--- + -- + 

k 

When we consider an even particle state, e.g., the vacuum, we take 

p(2k) = (+), p{2k + !) = {-) forfc^ioo (2.12) 

as the boundary condition. On the other hand, we take 

p(2k) = (T), p(2fc + 1) = (±) forfc^ioo (2.13) 

for an odd particle state. The ordinary approach of the Bethe Ansatz dismisses all 
the odd particle states, because the cyclic boundary condition on a lattice with even 
number of sites is discussed. Let us go back to the discussion of the q expansion. If the 
particle-number is greater than one, there is no a priori choice of the limit q = 0. The 
shift T 2 can fix only the total momentum (at q = 0), which is certainly not enough. 
So we need some new idea other than (i) and (ii) above. The idea is, of course, to 
use the quantum affine symmetries discussed in §1. What can one do if one has the 
symmetries of the Hamiltonian? The answer is to decompose the space of states into 
the irreducible pieces. This idea works very well if the irreducible decomposition is 
multiplicity free. As for the XXZ model in the anti-ferroelectric regime, this is not the 
case if we consider only the local symmetries and the non-local f/ 9 (sl(2))-symmetries. 

But if we consider the C/ g (s[(2))-symmetries this is what we actually get. In the case of 
the CTM Hamiltonian, the dccomposability of the physical space is rather trivial; the 
space V(A ) is already irreducible. This is one big reason why the CTM magic works 
so nicely. For the XXZ Hamiltonian, we should "decompose" V(A ) ® V(A )* a , (or 

equivalently EndQ( 9 )(V(Ao))) which is a level representation of ?7^(s[(2)), an object 
highly reducible and much more interesting than the irreducible modules. 



■)e lku 

■y ku 
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§3. Decomposition of Crystals 

At q = the Hamiltonian (4.1) is diagonal. Apart from the divergence of its eigenvalue, 
a vector of the form 

• • • <g) v pi 2) <2> f p (i) ® v p (o) ® «p(-i) <8» ■ • ■ 

for an arbitrary map p : Z — > {+, — } is an eigenvector of the Hamiltonian at q = 0. 
Among these, we consider only those which satisfy appropriate boundary conditions. 
We call them paths. The set of paths is equal to the product of two affine crystals, those 
corresponding to V(Ao,i) and V(A )* a . (See [2] for affine crystals). We also consider 
the q = limit of the creation operators tpj. (We borrow the commutation relation 
(7.106) from section 7.) We introduce the algebra generated by ip*, and show that it is 
isomorphic to the crystal of the paths. By using this isomorphism we give the q = 
limit of the n-particle eigenvectors. 

3.1 Crystal V 1 Recall that the highest weight module V(Ai) (i = 0, 1) has the crystal 
which are realized as paths [29] (see below). Let us denote it by V\ e $ t . Similarly, we 
have the crystal V l right for V(A l )* a . The crystal V° A of V(A ,i) <g> V(A )* a is the tensor 
product ~Pief t <8> bright m tne sense 01 crystals. As a set it is just a direct product 
V 1 = Vfeft x Vrighf equipped with the action of the modified Chevalley generators e,, fi 
according to the tensor product rule [25]. 

Practically the crystal V 1 is described as follows. An element of V 1 is represented 
as a path extending to both directions: 

P= (P( k ))kez> P(k) &{+,-}, (3.1) 
p(k) - (-) k for k < 0, p(k) = (-) k+l for k > 0. 

The labeling k is from right to left, k > 1 (resp. k < 0) for V t ^ t (resp. V^ ight ). 

Alternatively V 1 can be described as the set of arbitrary decreasing sequence of integers 
?!>•••>/„ with n = i mod 2. The two pictures are related via 

{k E Z | p(k) = p(k + 1)} = {h, • • • , l n }, h>--->l n . 

We call l\, ■ ■ ■ , l n 'domain walls' of the path p, and write . . . , l n ]] to represent p. We 
have the obvious decomposition 

V a = |J V(n), V 1 = |J V(n), (3.2) 

n:even n:odd 

where V(n) is the set of paths with n domain walls. 

There are two types of domain walls: adjacent (++) and ( ) pairs. By definition 

the spin s(p) of a path p is simply ((((++) — |( ))/2. The total weight of (3.1) has the 

form 

wtp = s(p)a\ — h(p)S. 
The h(p) £ Z is given in terms of the energy function H{e,e') 

H(+, -) = -1, H(e,e') = otherwise 
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as follows: 

Hp) = E k(H(p(k), P (k + 1)) - H(p l {k)Mk + 1))) • 

feez 

Here p £ P l and pi(k) = {-) k for k < 0, pi(fc) = (-) fc+ * for fc > 0. 
The rules for the action of , /j may be expressed as 

fo ■ (-) -> (+), /o : (-+) -> 0, 
A : (+) - H, A : (+-) -> 0. 

Their application is as follows (we shall describe fo). fo has no action on any ( — h) 

pair (a singlet) . Given a path p we first cut its left tail • • • — I 1 h and right 

tail — I 1 h • • • to make it a finite sequence. We reduce it further by the rule 

that each time there is an adjacent singlet pair ( — h) we drop that pair. For example, 

(••• — hH h H 1----) reduces to (+). The above procedure is not unique but the 

result does not depend on the particular choice of it, reflecting the coassociativity of the 
coproduct. If no spins are left then the path as an entity is a singlet with respect to fo 
and it is annihilated. Otherwise we necessarily have a sequence (H 1 — • ), and 

fel k 2 

it is mapped to (H 1 — • ). Namely, if fc 2 7^ 0, one and only one — is changed to 

fci+i fe 2 -i 

+ by the action of fo- Note that the configuration of the path (before the reduction of 

singlets), at the next left and the next right to this (— ) is (H ). This is changed into 

(+ H — ). Thus fo shifts a domain wall to the left and cannot cause two domain walls 
to coalesce. In the above situation eo produces (H 1 — • ). The rules for £1, f\ are 

fel-l + 1 

given by exchanging the roles of + and — . 

Consider first the case n = 1. The paths in V(l) have spin ±1/2 according as I 
is even or odd. They are a crystal of Aff(£?) according to Definition 2.2.3 of [2] with 
non-zero action (I £ 2Z) 

A:[[Z]]~[[/ + 1]], fo : [[I ~ 1]] ~ [[/]], le2Z, (3.3) 

and weights 

wt[[/]] - -(1/2)5 + (1/2H, wt[[l + 1]] = -(1/2)6 - (l/2)ai. 

Thus V(l) is a connected crystal. 

Now consider n — 2. Corresponding to (3.3) we have (with £ 2Z, l\ > I2) 

h ■ [[h - 1, h]] ~ [[h,h]] ~ WhM + 1]] (3.4a) 
fo ■ \[hM + 1]] ^ [[h + l,h + 1]] ^ [[h + l,h + 2]] ^ (3.46) 

and the weights of these triplets are respectively 

h + h 



2 

h + h + 2 



5 + ai, 0, —ai for (3.4a), 
5 + cto, 0, — a for (3.46). 
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So V(2) has an infinite number of disjoint connected components: each such component 
may be labeled by the fact that it is generated from the spin-1 path [[I, 0]] (I E 2Z> + 1) 
by (3.4). V(2) is only "half" of M(B) <g> M(B), the latter being labeled by a pair of 
integers without restriction. 

Generally, in V(n) the maximum spin of a path is n/2, and such a path has all (++) 
walls. Typical is 

[[h,l 2 , •■■,*„]], e2Z> + l (3.5) 

with wall labels alternating between even and odd integers and weight 

trn / 7ii ( l i + --- + l n+[i]) x , n 

Wt[[«i,i 2 ,---,'njj = ^ -<5+2 ai - ( 3 - 6 ) 

Here [n/2] denotes integer part. Because of the ordering of the wall labels, there arc in 
general only "1/n! times as many" connected components in V(n) as in Aff(B)®" . The 
precise formulation will be given in the next subsection. 

3.2 Paths as a quotient o/U£L M(B)® n Let Z be the Z-algebra generated by ip* (j e 
Z) satisfying the relations 

^+^ = 0, v; + i^ + ^+i^ = o, (3.7) 

for all j, k e Z such that j = fc mod 2. This algebra arises from the commutation relation 
(7.106) of the creation operators expanding then formally in q and z. 

ip* + (z) = J2r2 j z j + 0(q), 

<p*_(z) = J2r2 j+ iZ j + 0(q). 
jez 

Here we assume that the creation operators preserve the crystal lattice, i.e., there is no 
negative powers in q (see Appendix 4). We have also removed the fractional powers in z 
from tp±. 

Consider the Z- module ZV spanned by the set of paths V° U V 1 . We define an 
action of Z on TIP as follows. Take [[/Ji, . . . , l n ]] £ V(n). If n = we set 

m ]] = [[?•]]• 

If j > l\+n, then we define 

i>j[[h, ■ ■ -Jn}} = [[j - n,h, . . .,l n }]. 
H j = h + n or li + n — 1, we define 

v*[[/i,...,y] = o. 

Finally, if j < l\ + n — 1, then we define 

^[[Zi, . . . , /„]] - -< +n _i Vj[p2 ■ • ■ U] if J = ii + n - 1 mod 2, 
= -K+n^-i [fo ■ ■ ■ In}} if j # Zi + n - 1 mod 2. 
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From this action we see that the following is a linear base of Z. 

B = U™ =0 B{n), B(n) = {V& ■ ■ ■ ^ n ;ji - n + 1 > j 2 - n + 2 > ■ ■ ■ > j n }. 
There is a bijection uo from Z to TIP which maps B(n) to V(n). 

uj : B(n) -» V(n), (3.8) 
' " ' rfn ^ Ib'i - n + 1,J2 - n + 2, . . . , j„]]. 
This is obtained from the action of Z on the vector [[ ]] e "P(O). (Inparticular, = 

We make an affine crystal as follows. We define the weight by 

n 

wt(v<^--vi) = £wt(v>*j, 
fe=i 

w t(V>2.j) = -i* + wt (^2j+i) = -js - ^<*i- 

We define the crystal structure on B(l) by identifying it with AS(B) by ip*- = [[j]]. 
To define the crystal structure on B(n) (n > 1), let us consider the map 

U~ o Aff(B)® n - u~ (B(n) U -B(n)) U 0, 

Starting from the product tp^ ■ ■ ■ ip*^ with an arbitrary set of indices j = (ji, . . . ,j n ) G 
Z", we can modify it by using the rule (3.7) to ±^*, •••ipj, where j' k > j' k+1 + 1 for 

all k, or to 0. This process is compatible with the arrows in the crystal U^L Aff [B)® n 
in the following sense. Suppose that an index set j changes to j' by applying the rule 

(3.7) once. Suppose also that an arrow goes from tp^ ® • • • ® -0j„ to ^ ® • • • ® ipj , and 

i)% ®- ■ to ?A* ® • • • ® i>*, ■ Then, by case checking one can show that j changes to 

j' by applying (3.7) once. Therefore we can induce a crystal structure on B consistently 
from U~ =0 M(B)® n . 

The map lo from B{n) to V(n) is an isomorphism of crystals, i.e., it is a bijection 

and it commutes with /; and efj. This is not an isomorphism of affine crystals, because 
the affine weights in B(n) and V(n) differ by a mutiple of 5: 

wt(w(V£ ® • • • <S> ipjj) - wt(V'* 1 <S> • • • <2> tp* n ) = n 2 /A if n is even 

= (n 2 - l)/4 if n is odd. 

<?.<? TTie q = 2irm£ o/ i/ie eigenvectors 

We discussed in Section 2 that there is no a priori method to forsee what the q = 
limit of the eigenvectors of the Hamiltonian is. Here we give a prediction using the map ui 

(3.8) . Note that the definition (3.8) is based on the commutation relation of the creation 
operators which will be discussed in Section 7. 

Define 

r + (z) = £ z^* 2j , r- (z) = £ z^j+i 

jez jez 

We conjecture that the q = limit of the n-particle eigenstates (see Section 5 and 7) are 
given by 

wty e \(*i) •■■<,(*"))• 
The validity of this conjecture is checked for n = 1, 2 in Appendix 2. It is also consistent 
with the Bethe Ansatz calculation for n = 2h and s = h, h — 1 in Appendix 5. 
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§4. The vacuum state 

The vacuum state is the unique U' q {si{2)) -singlet in W under the boundary condition 
(2.12). To find the vacuum in W means to find a C/g(st(2))-linear map 

Q(q) -^V(Ao)®V(AoT a . (4.1) 
It is given by the canonical element; 

|vac)=5^v fc <8>v£. (4-2) 
k 

Here and v k are dual bases. Since V(A ) is infinite dimensional, the above sum is 
an infinite sum. Therefore we need a completion of V(Aq) ® V(Ao)* a in order that 
the vacuum vector actually belong to it. For our purpose the way of completion is 
not very important. We take the largest one, i.e., the direct product of all the spaces 
V(A )\ <S> F(A )* a where A and \i runs the weights of V'(Ao) and V(A )* a respectively. 
We do not use any particular notation for the completion. We just use V(Ao) ® V(Ao)* a 
as being completed. It is convenient to use another realization of V^(Ao) ® V(Ao)* a . 
Namely we use the canonical isomorphism 

V(A ) ® V(A y a ~ End Q(9) (V(A )) . 

(Strictly speaking, the right hand side is the set of Q(q)-linear homomorphism from 

V^(Ao) to the completion of V(Aq).) The f/ ? (sl(2)) -action on the left hand side is given 
by the coproduct. Then the action on the right hand side is given by the adjoint action 
which we denote by ad x; suppose that 

/ e End Q(9) (V(A )), v e V(A Q ), x e U q (sl{2)) 

and 

A(a:)=X;4 1) ®4 2) - 
k 

Then we have 

((ad x)f)(v) = J2x^f{a(x^)v). 
k 

If / = id, then (ad x)f = e(x) (e : the counit). Here we used the special choice of the 
left module structure V* a in (4.1). Therefore id £ Endq( g ) (V(Aq)) generates a singlet. 
Or equivalcntly, the canonical element (4.2) actually realizes the vacuum. In Appendix 
4 we give some explicit computation of the embedding of the canonical element in W. 

One may raise a question about the translational invariance of the vacuum. Because 
of the splitting of the whole line into the right and the left pieces, it is not obvious that 
the definition of the vacuum is independent of the choice of the position of the splitting. 
Let us distinguish the vacuum in V(A ) ® V(A a )* a and the vacuum in V(Ai) ® V'(Ai)* 
by denoting the former by \vac)o and the latter by \vac)\. We want to prove 



T\vac)i = \vac)i-i i = 0,1. 



(4.3) 
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In order to prove this we need the interpretation of the translation T in our mathematical 
picture. For this purpose we use the vertex operators 

$ : V(Ai) — » y(Ai_i) ® V 

and 

** : V <g> V(A 4 )* a — ► F(A 1 _,)* a . 
We introduced $ in §2 with the normalization 

$(w Ao ) = «Ai ® (-) + •■• , *(«Ai) = «Ao ® (+) + ••• • 

We define 'J* as the intertwiner with the normalization **((—) <g> u Ao ) = u% i + ■ ■ ■ and 

**((+) <8> u Al ) = u Ao H , where w Aj is the lowest weight vector of V(Ai)* a dual to u Az . 

We define 

T : V(Ai) ® t/(A 4 )* a — V(Ai_i) ® F(Ai_i)* a 

by the composition 

V(Ai) ® t/(A 4 )* Q *^ d F(A 1 _ l ) ® y ® ^(AO^^VlALi) ® t/(A!_ t )* a 

up to a constant multiple. This constant is determined in §7 and the proof of the assertion 

(4.3) is also given (see (7.4), Proposition 7.1). Once we establish the statement (4.3) then 

Hxxz\vac) t = (4.4) 
follows immediately, because by the definition it is obvious that 

d\vac) l = 0. (4.5) 

Before ending this section, we wish to comment on [31], where it was first conjec- 
tured, on the basis of direct perturbative calculations up to a low order in q, that the 
vacuum vector of the XXZ Hamiltonian is also an eigenvector of (1.4). In our algebraic 
scheme, this follows directly from (4.5). Our scheme also makes it clear that the excited 
states of the XXZ Hamiltonian are not eigenstates of d. 

Before ending this section, we will make a comment on [31]. In this paper a conjec- 
ture is given that the vacuum vector for the XXZ Hamiltoian is also an eigenvector of 

(1.4) . This is trivially correct in our scheme by (4.5). On the other hand, it was rather 
surprising when it was first pointed out in [31]. We note that the excited states of the 
XXZ Hamitonians are not eigenstates of d. 
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§5. Particle Picture 

By a particle we mean a finite-dimensional L^(sl(2))-module consisting of eigenvectors 
of the XYZ-Hamiltonian. ^ 

Starting from the infinite tensor product W, we have argued how the U q (sl(2))- 
module 

V(A ) ® V(A y a ~ End Q(g) (V(A )) 

is embedded therein, and how the vacuum (i.e., the zero particle state) is understood 
as a vector of the latter. Now, guided by the crystal decomposition of the even and the 
odd paths, we reach the following; Modulo statistics to be discussed later, the particle 
picture for the AXZ-Hamiltonian in the anti-ferroelectric regime is 



T = 



V, ® • • • ® V Zl du n ■ ■ ■ dux 



(5.1) 



symm 



where Ui is a quasi-momentum of the 2-dimensional C/^(sI(2)) -module V Zk [z^ — e lUk ). 
We call T the Fock space. (For the meaning of the symmetrization, see (7.6).) The 
even particle (n: even) are contained in the even sector V(Ao) ® ^(A )* a , and the odd 
particles (n: odd) in the odd sector V(Ai) ® V(Ao)* a . To see this we want to find a 

U q (sl(2))-lmeax map 

V Zn ® • • • ® V Zl - V(A or i) ® V(A y a , 
or equivalently, a f7^(s[(2)) -linear map 

V Zn ® • • • <8> V Z1 ® V(Ao) -» V^(A or i). 

Therefore, the problem reduces to finding the vertex operators 

: ® V(Ai) - F(Ai-i) (i = 0,l). 

The existence and the uniqueness of such vertex operators are given in [D JO] in a general 
setting. In Appendix 2, we give some explicit calculation of the above vertex operator 
in terms of the global base of Kashiwara. 

Let us argue the physical content of our particle picture. The AXZ-Hamiltonian 
possesses the infinite hierarchy of the abelian higher order Hamiltonians, and the infi- 
nite dimensional non-abelian symmetries of C/g(s[(2)). These are the symmetries of the 
AAZ-Hamiltonian in the strict sense, i.e., they commute with Hxxz- The abelian sym- 
metries do not change the eigenvectors, and the non-abelian symmetries do not change 
the eigenvalues. The Lorentz boost e £d (e : a scalar parameter) actually changes the 
energy, but never creates new particles nor annihilates them. Now we introduce the 
third symmetry of the AAZ-Hamiltonian, the dynamical symmetries which create and 
annihilate particles. 

Define the creation operator 

<p* ± {z) : V(Ai) ® V(A y a -> ViA^i) ® V(A )* a 

by 



<p±(z)(v ® v*) = ®(z)(v± ® v) ® V*. 
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Then, it is easy to see that <p±{z) acting on an n-particle state create an (n + l)-particle 
state. In this way, the Fock space T is embedded in (V(A ) © V(Ai)) <g) V(A )* a . 

In §6 and §7, we give a mathematical treatment of the creation and the annihilation 
operators. Here, we give rather heuristic discussions on several points. 

We argued the embedding of V(A ) ® y(A )* a to the infinite tensor product. The 
vectors G(p) in V(Ao) are expanded in terms of the paths \p), where we consider the 
latter as vectors in the infinite tensor product. A question arises: Is it possible to expand 
the paths in terms of the vectors in V(A ). The answer is no. The infinite matrix which 
expresses the trasition from \p) to G(p) is of the form 1 + qA\ + q 2 A 2 + ■ ■ ■■ If we invert 
this, we get formally 1 — qAi + q 2 (A\ — A 2 ) + ■ ■ ■■ But A\ has divergence on the diagonal. 
So, the transition matrix is not invertible. This means that our definition of the creation 
operators does not apply to the paths. Suppose we were to apply the creation operator 

to the bare vacuum (• • H 1 •••). Since the action of the creation operators changes 

only the left half, the right half is unchanged. This contradicts with the fact that the 
creation operator satisfies the proper commutation relations with the shift operator (see 
section 7). 

We conjecture that the creation operators preserve the crystal structure in the fol- 
lowing sense. The vacuum vector embedded in the infinite tensor product is expanded 
in power series of q. The conjecture is that the particle states created by the creation 
operators upon the vacuum also have the same property. This is remarkable because 
the vertex operators used in the definition of the creation oeprators do not preserve the 
crystal lattices. In fact, if they do we have the following contradiction. At q = the 
vacuum considered in V(A ) (g> y(A )* a reduces to ma <8> u\ g , where u\ is the highest 
weight vector in V(Ao) and u* A0 the dual lowest weight vector in V(A )* a (see Appendix 
4). So, if the crystal lattice were preserved by the vertex operators, then to get the q = 
limit of a one-particle state, we only have to apply the vertex operator to wao- But, this 
breaks the translational covariance which is expected for the one-particle state. In fact, 
we will see in Appendix 4, that we have contributions to the q°-teim in the one-particle 
state from the higher order terms in the vacuum state. 

Now we come to a more subtle point. In Appendix 4, we computed 1 and 2 particle 
states by applying the vertex operators successively. In the computation of 

$(z 2 )(w_ ® $(zi)(w+ ®u Ao )) 

(which is a part of the computation of ip*_{z 2 )Lp* + {zi)\vac)) , we find terms having poles at 
q = 0. But these terms seem to be summed up to a meromorphic function in zi, z 2 and 
q that has no pole (actually has a zero) at q = 0. So, even though the vertex operator 
<f>(z) : V z ® V(Ai) — ► V(Ai-i) does not preserve the crystal lattice, the particles created 
by <p±(z) may not (and, in fact, do not) have poles at q = 0. Namely, the crystal picture 
survives. The necessity of this summation also tells that the Fourier components of the 
creation and the annihilation operators are not equal to those of the vertex operators. 
We will further discuss this in Section 7. 

Let us mention a few words about the statistics of our particles. The tensor products 
V Z2 ® V Zl and V Zl ® V Z2 are different but isomorphic. They are intertwined by the R- 
matrix. In fact, the embedded images of these tensor products in V(Ao) ® V(Ao)* a are 
equal. This follows from the uniqueness (up to scalar multiple) of the vertex operator 
V Zl <g> V Z2 <S> V(A ) — > V(A ). So the key point is to determine this scalar multiple. We 
will give the answer to this question in §6 by using the quantum K-Z equation. 

Do particles with spin higher than ^ exist? Let (j G |Z) be the (2j + 1)- 

dimensional U q (sl(2)) -module, and V Z ^ its J7g(s[(2))-cxtension. If we understand this 
question as the existence of the vertex operator 



V® ® V(A ) -» V(A ) or K(Ai), 
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one can show the non-existence by an argument of crystals. If we treat the problem 
honestly, we should start from the finite-lattice Bethe Ansatz and examine the thermo- 
dynamic limit. There are papers [11] for the XXX-case (A = —1) and [9] for the 
XXZ-case (A < —1), which are in support of our picture. 

Let us consider the dual space V* of V. The AXZ-Hamiltonian formally acts on 
the infinite tensor product of V*, too. So, we can formulate the theory in a completely 
analogous way by using right modules. In that case, the (even) particle picture will be 

developed on the tensor product of two right modules V(Ao)* <g> I^(Ao)** a . Let us, in 
general, define a natural pairing ( ) 

(V(A')* ® V(Ay* a ^) ® (V(A') ® V(A)* a ) Q(q), 

where A and A' are arbitrary dominant integral weights. The left f7 g (s[(2))-module 
V(A') ® V(A)* a and the right U q (sl(2)) -module V(A')* ® V{\)** a have the canonical 
pairing ( , ) induced from the pairing between V(A') and V(A')* and the pairing 
between V(A)* and V(A)**. This pairing ( , } satisfies 

{fx, 9) = {f,xg}, 

for / G V(A')* ® V(X)** a , g G V(A') <g> F(A)* a and x G ?7 g (st(2)) . Note that there is an 
isomorphism of J7 9 (s[(2))-modules V(A)** a ~ ~ y(A)** a given by 

(v(A)** a ~ 1 ) x 3»h g - 4 ^> A ) w g (y(A)** a ) A . 

We define ( | ) by setting 

(u' <8>v/\u®w) = ( q - 4 (P^(w')) u > ® W < 7U(E)W ^ 

for weight vectors u',w',u,v. If we regard / G V(A')* ® V(A)** a 1 (rcsp. <? G V(A') ® 
F(A)* a ) as an element of Hom(V(A), V(A')) (rcsp. Hom(V(A'), V(A))) then we have 

{.f\g)=tr v{A) (q-^fog). 

Obviously this pairing satisfies 

(fx\g) = (f\xg) forxeC/ 9 (s[(2)). 

In particular, we have 

(vac\vac) = tr y(Ao) (<T 4p ) , 

which is the specialized character for V(Aq). 

Getting the character expression for the one-point function of the six-vertex model 
[2,23] is a simple corollary of this formula. Define a non-local operator 

j>k 
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on V (Ao)®V(A )* a . We understand the meaning of the infinite sum by the normalization 
To(|Ao) £g> | — A }) = 0. The one-point function 

(vac\Proi(To=rn)\vac) 
{vac\vac) 

of the six-vertex model is by definition the expectation value of the projection operator 
to the t = m cigenspace. Therefore, it is given by 

P(m) = try(A ° } -C * 

where V(A ) m is the spin m subspace of V(Ao). 

By a similar construction, the dual Fock space T* is embedded in (V(Ao)* © 
^(Ai)*) ® V(A )** a_1 . The bilinear coupling between V(Ai) <g> V(A )* a and V(Aj)* ® 

V(A )** a ~ induces a non-degenerate coupling between T and T* . The n-particlc states 
V Zn <g> • • • <S> V Zl in J 7 and the m-particle states V*, <S> • ■ • <g> V*,^ in T* are orthogonal unless 

n = m and {zi} = {z^}. 

Finally, we will consider the annihilation operator. As we defined the creation oper- 
ator in the frame of left modules, we define the operator which create finite-dimensional 

right U' q {s{{2)) -modules, i.e., the dual Fock space, in V(A )* ® V(A )** a . It is given 
by means of the vertex operator of the form 

Let v J <G Vj* be the dual elements to v± ; (u* ,v E >) — 5 es > . Define 

<p ± {z) : V{A^* ® V(Ao)** a " -> V^(Ai_i)* ® y(A )** a " 

by 

</?±(z)(w* <g> w) = $*(z)(w± <g) w*) (g) w, 

where w* G l^(Aj)* and u; G V(Ao). We define the annihilation operator ip±(z) acting 
on T by the dual action of tp±(z) : T* — > T* . Since </?±(z) creates (n + l)-particle states 
from n-particle states in T* , y±(z) annihilates (n+ l)-particle states to n-particle states 
in T . 
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§6. Vertex operators 

So far we have argued how one can treat the elementary excitations in the anti- 
ferroelectric XXZ model in the framework of representation theory. The motivations 
being given, we now turn to the mathematics of vertex operators. 

6.1 Notations Let us restart by fixing the notations. Thus let P = ZAq © ZAi © Z5 

and P* = Zho © Zhi © Zd be the weight lattice of st(2) and its dual lattice respectively. 
We have (A it hj) = 6 ijt (A*, d) = 0, (S, hi) = and (S, d) = 1. We set a x = 2A X - 2A , 
a Q = 5 — ai and p = A + A x . We normalize the invariant symmetric bilinear form ( , ) 
on P by (ctj, Qj) = 1. Explicitly we have (A,, Aj) = SnSji/A, (Aj, 5) = 1/2 and (5, S) = 0. 
We shall regard P* as a subset of P via ( , ), so that 2oti = hi and 4p = /ii + 2(i. 

The quantized affine algebra U — U q (s{{2) ) is defined on generators e,, /j (i = 0, 1), 
q' 1 (h e P*) over the base field Q(<j). The defining relations are as given in [2,25], e.g. 
[ei,fj] = SijiU - ^/(q - q- 1 ) where U = q h \ We let U' = U' q {s[{2)) denote the 
subalgebra generated by ej, fi and ti (i = 0, 1). We shall take the coproduct A to be 

A(e 4 ) = e 4 © 1 + t t © e», A(/i) = / f © t^ 1 + 1 © ft, 
A(q h ) = q h © q h (h E P*). 

(6.1) 

Accordingly the formula for the antipode a reads 

o(e») = -*r le ^ = a (<?' 1 ) = 'T' 1 - 

6.2 Modules Given a left [/-module M we write its weight space as M v = {v G A'/ | 

q h v = q^v yh e P*}. For u e M u we write wt (u) = v. Suppose M = ® V M V , and 
let 4> be an anti-automorphism of U. Then the restricted dual M* = ® V M* is endowed 
with a left module structure M*^ via 

(u, xv) = (<j>(x)u, v) for x E U, u e M, v G M* . 

We have M ~ (M*^)*^ 1 (canonically). Similar convention is used for right modules 
and [/'-modules. Taking (f> — a we have the canonical identification 

Homy- (L, M © N) = Homy- (M* a © L, N) 
Homy/ (L ®N,M) = Hom^ (L, M © A^* a ) 

(6.2) 

for left modules (for right modules we replace a by a -1 ). 
In what follows we set 

A = Ai, fi = Ai_j for i = or 1. (6.3) 

As before let V(X) (resp. V r (Xj) denote the integrable left (resp. right) highest weight 
module with highest weight A. To distinguish the left and right structures we shall use 
the bra-ket notation (u\, \u) for vectors in V r (X) and V(X), respectively. We fix nonzero 
highest weight vectors {u\\ G V r {\), \u\) G V(A) once for all. Then there is a unique 
non-degenerate symmetric bilinear pairing V r (X) x V(X) — > Q(q) such that 

(ua|«a) = 1, (uiclu') = (u\xu') for any (u| G V" r (A), |u') G V(A). 
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We shall also consider the simplest two-dimensional [/'-module V = Q(q)v + + 
Q(q)v_ given by 

eiv + = 0, e\V- = v+, fiv+ = V-, fov- = 0, t\v± = q ±1 v±, 
eo = fi, fa = ei, t = t^ 1 on V. 

For fixed meZwe equip V <g> Q(q)[z, z" 1 ] with a [/-module structure by letting 

d act as e; ® (zg™)' 5 * , /; act as /; ® (zq m )~ Sm , 
the weight of v± (g> z n = n5 ± (Ai - A ). 

The resulting [/ -module will be denoted by V zq m . (This conflicts with the notation for 
weight spaces, but the meaning will be clear from the context.) 

Let v± be the basis of V* dual to v±: (vi,v*) = 5ij. One checks readily that the 
following give isomorphisms ('charge conjugation') of [/-modules 

C± : V zqT 2 -^V z * a±1 (6.4) 



' zq 

C±v + = v*_, C±v_ = —q ±x v* 



+ • 



6.3 Basic vertex operators By vertex operators (VOs) we will mean the intertwiners 
of [/'-modules of the type 

$f : V(X) — >V(jt)®V, (6.5) 
: V{\) — >V®V(ji). (6.6) 



Here V(fJ,) — \\ v V(p) v is a completion of V(fj,). Since the weight is preserved modulo 6, 
for given v € ^(A)^ one can write v = X)ngz( u +,« ® v + + u -,n ® u_) where u±, n G 
^ / (a i )i>=f(Ai-a )+"<5- (Because the weights of V(/x) are bounded from above, u± in = for 
n large enough.) Thus one can define the weight components ($^ V ) ± , ($[ M ) ± by 

*f = e (*f) ± ,„®^. *r"= e ^®(*r") ±lB . 

n£Z,± n£Z,± 

(*f) ± , n , : f(A), — V(^ {Al _ Ao)+nS . (6.7) 

We shall fix the normalization as follows: 

$5t V |«A> = M® v =F + -". ( 6 - 8fl ) 

$^|«A>=%®|«m> + -"- ( 6 - 85 ) 

Here t>_ (resp. t> + ) is chosen for A = Ao (resp. A = Ai). In (6.8a) • • • means terms of the 
form \u) <giv with \u) £ V(fi)^, and similarly for (6.86). The existence and uniqueness of 
such VOs are shown in [1,4]. 

The VOs can be equivalently formulated as intertwiners of [/-modules of the form 

[1] 

*?(z) = $f(z)A- A \ = i^wA-^, 
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Here we set A a = (A, A + 2p) /(k + h v ) where k = 1 is the level and ft, v = 2 is the dual 
Coxeter number; explicitly Aa = 0, =1/4. The right hand sides of (6.9a), (6.96) 

mean e.g. V{p)®V z = 0„ Yl^ V(fi)^ ® (V z ) v -£. The fractional powers are so designed as 
to put the q-KZ equation in neater form, see below. 

By abuse of notation we let d € End(V(A)) denote the operator 

d\u) = (d,u)\u) |u)eV(A)„. (6.10) 
It is easy to see that [d, (^^ V ) ± J = n (^\ V )± „ an d hence that 

(d ® id)4>f (z) *f(z)d =-( z ±-A, + A A )$f (z). 

Similar relation holds for $^(2). 

Remark. In [2,4] the coproduct of U is chosen to be 

A_(e 4 ) = e 4 ® tr 1 + 1 <g> e i; A_(/j) = /, <g> 1 + U ® 
A_(q' 1 ) =q' l ®o' 1 (heP*). 

The present formulation using the coproduct A = A + (6.1) is related to the references 
above as follows. 

Let M, N be U- modules such that wt (M) C A + J2 Za n wt ( N ) C Mo + E z "i for 
some Ao,/io G -P. We define operators /3m, 7mat by 

V = g 2(A,/i)-2(Ao,M) u0 ^ u g g 

Then A+(x) = 7mjv° A_ (x)o^ n (x e U) and/3 M ®/3Af = Pm®n°1mn = 1mn°Pm®n- 
We extend 7m jv also to M®N . It is known [25] that 

(i) (L, £?) is a lower crystal base of M if and only if ((3m(L), 0m(B)) is an upper crystal 
base of M. 

Suppose Mi have lower crystal bases (L 1 ™ , B\ ow ) {i = 1, 2, 3), and set L" p = /3 M ,(i' ow ), 
= /3m s (B] ow ). For $ low : M x -> M 2 ®M 3 we put $ u p = 7m 2 m 3 °$ 1ow Then we have 

(ii) $ low a; = A_(a;)$ low if and only if $ up x = A + (a;)$ up (xeU), 

(hi) $ up /3 Ml = /3m 2 ®/3m 3 $ 1ow . Hence $ 1ow (l5l >w ) C Z, 2 ow §£ 3 ow if and only if $ u p(L" p ) C 

Similar statements are valid for the intertwiners of the type * : Mi ® M 2 — > M3. 

6.4 Variants of vertex operators The identification (6.2) along with the isomorphisms 
(6.4) gives rise to the following variants of intertwiners. 

Type I: 

$f : V(X) —> V{p) ® V, (6.11a) 

$ A y :V(p)®V — V(A), (6.116) 

*Vm : V" <8> V(^)* a — » V(A)* a , (6.11c) 

: V(A)* a — ► V <8> V(^)* a . (6.11d) 
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Type II: 

&X»:V(\) — >V®V(fi), (6.12a) 

fyr-.VQVOjL) — ?(A), (6.126) 

^ :V(fi)* a (g)V — >V(X)* a , (6.12c) 

: V(\)* a — > V(n)* a Cg) V. (6.12d) 

That is, (6.116) is obtained from 

— ►^(A)®V*°~\ 

and (6.11c), (6.11d) are transpose of 

V(A) — > V(fj,) <S> V"(/x) <8> F* a_1 — >T>(A) 

respectively. The case of Type II is similar. We define 

^ xv (z)(u ® „) = (id ® (v, ))$f "(*)«, $f *>) = (id ® C+)$f (zq- 2 ), 

(6.13a) 

^ vx (z)(v ® «) = ) ® id)$r a_v (^K = (c- ® id)$^^g 2 ). 

(6.136) 

This implies the normalization 

#£ v (*)(I«a> ® v±) = T« ±1/4 g 1/2 |«„> + • ■ ■ , 

$^ A (z)( V± ® | UA » = T« ±1/4 3- 1/2 |«„> + ' • •> 

where the upper (resp. lower) sign is chosen for A = Ao (resp. A = Ai). For example 
Z -A„+A x $i* v ( z j is a [/(g)Q(q)[ ZjZ -i]_li n ear map 

V(\)®V z —>V(fi) ® Q(q)[z, z- 1 ] 

where on the right hand side x & U' acts as x®id and q d acts as q d <g>q zd / dz , (q zd / dz f)(z) — 

We have distinguished the two types of intertwiners. The main difference is that the 
type I operators preserve the crystal lattice (see Sect. 6.7), while the type II operators 
do not. We shall see explicit examples of the latter phenomenom in Appendix 4. 

6.5 Two point functions Frcnkel and Rcshctikhin [1] showed that the correlation func- 
tions of the vertex operators satisfy a q analog of the Knizhnik-Zamolodchikov (q-KZ) 
equation. For our subsequent discussions we need mostly the case of two point functions 
for various combinations of VOs. The general case will be discussed in Sect. 6.8. 
Let ^(zi,z 2 ) be one of the following correlations: 

(u x \^(z 2 )^( Zl )\u x ), (6.14a) 

(^<&^ 2 )<I>f^i)K>. (6.146) 

(u x \^(z 2 )^( Zl )\u x ), (6.14c) 

(u x \^ x (z 2 )^( Zl )\u x ). (6.14a-) 
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Clearly 

*(zi,z 2 ) e V ® 1/ ® (z 1 /z 2 ) A -- A -Q( <z )[[z 1 /z 2 ]]. (6.15) 

The subscripts of V in (6.14a)-(6.14d) indicate the tensor components. In (6.15) the left 
V is Vi and the right one is V 2 . Thus, for example, if we replace in (6.14a) by fi\u\) 
then the result becomes 

(ft ® 1 + tr 1 ® /,) (uaI^^)*^ 1 ^!)!^) e V <» V. 

(Note that in the discussion of the embedding of V(A) into V®°° in Sect. 3 and in Sect. 
6.8 we use the opposite ordering of the components.) In what follows we introduce the 

element 

q hi/4 i n 

U and extend the base field Q(q) by adding g 1 / 4 . (We could avoid using 
fractional powers, but the formulas would become slightly more cumbersome.) 
We need also prepare the R matrix. Define R(z), R + (z) by 



R(z)v± <S>v± = v± ® v±, 

R(z)v + ®v- = 7T-qv + ® w_ + r z))- ® v+, 

1 - q z z 1 - q z z 

1-q 2 l-z 

1 — q z z 1 — g z z 



R(z)v^ ® v+ = — v + ® i>_ + — qv ~ ® v +- 



(6.16) 



i?+(z)=p(z)i?(z), p(z) = a- 1 / 2 - '''^ 



(z) 00 (g 4 z) c 



where we put 

oo 

(z)^ = (z;g 4 )oo, (zjp)^ = ]J(1 - zp J ). 

j=o 

Let further P £ End(F®F) be = u'®u. Then PR(z 1 /z 2 ) :V Zl ®V Z2 ® 

is an intertwiner of [7- modules, and i? + (21/22) is the image of the universal i? matrix of 
U in End(T4i ® Vz 2 )- The scalar factor p(z) is determined by the argument in [1]. 
With these notations the g-KZ equations read as follows: 

*(g 6 z l7 z 2 ) = A(z 1 ,z 2 )V(z 1 ,z 2 ), *(g 6 z 1: g 6 z 2 ) = {q~* ® q~*) *(z u z 2 ), 

(6.17) 

where (j> = 4A + hi, Ao = 0, Ai = /ii/4, and 

A(zx, 2 2 ) = R + (q 6 Zi/z 2 )(q-' t> ® l) for (6.14a), 

= (g~ 2 ^ ® l)i? + (g 5 z 1 /z 2 )(g-^ +2X ® l) for (6.146), 
= (g 2X ^ ® l)i?+(gz 1 /z 2 )(g- 2X ® 1) for (6.14c), 
= (g- ® l)i? + (zi/z 2 ) for (6.14d). 

In the present case the property (6.15) and the normalization (6.8a)-(6.86) specify 
the solutions of (6.17) uniquely. We list the answers below. 
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A= An 



(z 1 /z 2 )- 1 / 4 vl/(z 1 , z 2 ) = (f Zl/z2 \°° (v- ®v + - qv + ® «_) for (6.14a), 
= ( q J* 1 /* 2 ] 00 ( v _ ®v+- qv+ ® v_) for (6.146), 



(g 3 Zi/z 2 ) c 

(gzi/z 2 ) 00 
(g-^i/^c 

(^l/z 2 )oo 
(g- 2 Zi/z 2 ) c 



-(w_(giw + — g 1 v + ®i;_) for (6.14c), 
-(w_ ® w + — g _1 w+ ® u_) for (6.14d). 



A = Ai: 



(zi/z 2 ) 1/4 *(zi, z 2 ) = i 9 / 1 ^ 2 ! 00 (t>+ ® «- - gzi/z 2 v_ ® w+) for (6.14a), 

(g zi/z 2 j 00 

- (« ^/«2)oo ^ + ^ y _ _ qZl / Z2V ^ w+) for ( 6 .i45) 7 



(g 3 zi/z 2 ) D 

(gfl/fgjc 

(g _1 zi/z 2 ) 



-(«+ ® v_ — g 1 zi/z 2 w_ (g) u + ) for (6.14c), 



oo 

(Zl/Z2)oc -(«+ ® v_ - q- 1 z 1 /z 2 v- ® v+) for (6.14d). 



(g 2 zi/z 2 ) c 

6.6 Commutation relations The general theory of g-diffcrence equations tells [32,33] 
that the n point functions of VOs can be continued meromorphically to the entire space 
(C x )™, apart from overall powers or logarithms in Zj. In our case this is apparent from 
the explicit formulas. It follows that the same is true of all the matrix elements of 
compositions of VOs (see the proof of Proposition 6.1 below.) 

From the knowledge of the two point functions one can derive the commutation 
relations of VOs [1]. To write down the relations which will be used later, we need to 
modify the scalar multiple of the R matrix and define 

R vv (z)=r (z)R{z), R v *v>(z) = (C_ ® C-)R vv (z)(C- ® C_)~\ 
R vv ,(z) = (id® C_)i? yy (zg- 2 )(id® C_)~\ 



(6.18) 



Here 



-i/2 r M (z-'utfzu ^Ah+^MzA ) 



° ( ) " WooW*- 1 )- - (6 ' 19) 

with z = q- 2 P/ l7T and r p (a;) = (p;p)oo/ (f> x ;p)oo(l — p) 1 ~ x denoting the g-gamma function. 
We have the unitarity and crossing symmetry: 

R vv {z)PR vv {z- 1 )P = id, 
(Rwiz)- 1 )* 2 = -R vv *(z). 
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Notice that Ryy{z), Ryy (z) and Ryv (z) have no poles in the neighborhood of \z\ = 1. 

In the following theorem, we list the commutation relations of VOs of type (6.12a)- 
(6.126). The commutation relations and the holomorphy properties are in the sense of 
matrix elements. For example (6.20c) below states that for each (u\ G V r (A), \u') G V(X) 
the following hold as meromorphic functions in zi, Z2 times {z\j z^^ 1 / 4 : 

(u\^ x (z 1 )^ a ~ 1 "(z 2 )\u') = R vv ,{z 1 /z 2 ){u\^ 1 \z 2 )^(z 1 )\u'). (6.20s) 

Proposition 6.1. 

(i) The following commutation relation holds: 

^\ Zl )^{z 2 ) = Rv^/z^l^z^^z,), (6.20a) 
^ 1 \z 1 )^^{z 2 ) = Rv-vMI^T ^i^T'^izi), (6.206) 
^(^i)*?" 1 ^) = Rv 1 v 2 ^z 1 /z 2 )^ 1 \z 2 )^(z 1 ). (6.20c) 

(ii) In the neighborhood of\zi/z 2 \ — 1 both sides of (6.20a)-(6.206) are holomorphic, 
while (6.20c) has a simple pole at z\ — z 2 . The residue of the latter is given by 

Rcs z=l ^ x { Zl )^Y ^{z 2 )dz = gxid v{x) ® (v+ ® v* + + w_ ® v*_ ) , (6.20d) 
where z = Zi /z 2 and 

with the sign — (resp. +) being fcaicen for A = A (rcsp. A\). 

Proof. The argument being similar, we concentrate on (6.20c). 

If (u\ = (u\\ and \u') = \u\) the assertions (i),(ii) follow from the explicit formulas 
(6.14) and (6.13a)-(6.13b). Suppose \u') = \xu") with x G U'. Then the intertwining 
property of vertex operators entails 

P{u\^ x { Zl )^ a ~^{z 2 )\xu'') 

= J2x(i)®X(2)P(ux {3) \^\ Zl )^ a ^(z 2 )\u") G VT" 1 ® VI. (6.22) 

Here A^ 2 ^(a;) = ^^(l) ® ^(2) <8) X( 3 ) . Analogous formula holds for {u\ — (u"x\. Since the 
action of U' on V* 2 a ® V Zl involves only powers of Zi, the analyticity property follows 
by induction on the weight of (u\, |it'}. Because PR(z) is an intertwiner the relation 
(6.20c) is unchanged in the process. To see (ii) note that if we take the residue of (6.22) 
then the left hand side reduces to 

= P{ux\^\z 1 )>S> V f~ 1 »{z 2 )\u"). 

This is because the vector w = v*L ® v + + t>* ® V- G V* a ® V z belongs to the trivial 
representation. We then obtain by the induction hypothesis (ux\u")g\w = {u\xu")g\w 
as desired. This completes the proof. □ 
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We shall also need the following commutation relations, which can be proved in a 
similar manner. 

R Vl v 2 (zi/z2)^(ziM V2 (z2) = -^(z 2 )^ Vl (z 1 ), (6.23) 
< lA (^)*fW - < 2 (^r^i) x ((^ 2 ) 1/2 )|[^fj, (6-24) 

$r" lA (, 1 )< 2 (^) = ^(^)$F" V (^) x (-(^iA 2 )- 1/2 )||g^. (6.25) 

Here we set 

e(z) = {zutfz-'utfu. 

It is known [4] that the composition 3>^y ° $^ is convergent in the g-adic topology 
and is proportional to the identity. The proportionality scalar can be determined from 
the two point function by setting z\ — z 2 . Using the results (6.14) we find 

*£v(*)o$$; v (*) = -i-xid (6.26) 

where g\ is given in (6.21). 

6.7 Crystal lattice Recall that on an integrable [/-module one has the operators e" p , 
/" p , e l ° w , fl° w and the notion of the upper/lower crystal lattice. [25]. The finite dimen- 
sional module V has L = Av + ® Av^ as a (both upper and lower) crystal lattice. 

The integrable highest weight module V(X) has the standard lower/upper crystal 
lattice L low (X), L up (\) given as follows [25]. Let <p denote the anti-automorphism of U 
given by 

<p(ei) = fi, <p(fi) = ei, p(q h )=q h , (6.27) 
and let ( , ) be the unique symmetric bilinear form on V(A) such that 

(|ua),|«a» = 1, {\xu),\u')) = (\u),\ip(x)u')). (6.28) 

Then 

L low (X) = the smallest A-module containing \u\) and stable under f\ ow , 
L up {\) = {ue V(X) | («, L low {\)) c A}. (6.29) 
L low (X) = {ue V(X) | (u, L u p(X)) C A}. (6.30) 

Here A = {/ G Q(<?) | / has no pole at g = 0}. We have further 

L up {X) v =q {x - X) - {v - l,) L low {X) l/ . 

The crystal lattices of integrable lowest weight modules are defined similarly by 
replacing \u\) by the lowest weight vector and f° w by e l ° w . Let T\ e End(V(A)) denote 
the linear map 

T x u = q( 2p ' x - u) u for ti e V(A)„. 



Diagonalization of XXZ Hamiltonian 31 

Proposition 6.2. The upper crystal lattice L* aup (X) ofV(X)* a is characterized as 

L* aup {X) = {ue V(X)* a | (u,T x L up (X)) C A}, (6.31) 
where { , ) denotes the canonical pairing. 

Proof. Let r( denote the map V(X) — > V(X)* a given by (rj'(u),v) = (u,v), and set 

»?(u) = {-l) ht ^- v )q("'"+ 2 P)-(^+ 2 P) ri '{u) u e V(A)„. (6.32) 
Here for £ = m a + m\Ci\ we set /it(£) = m + mi. One can verify that 

»?(eiti) = fiV{u), V(fiu) = eiT}(u), n(q h u) = q~ h n(u). 
This implies f](fl° w v) = £ l ° w n(v), f](e l P w v) = f\ ow ri(v), and hence 

L* alow (X) =j 1 (L low {\)). 
Using the relation (L* aup (X)) _ v = q^*) -<-">") (L* alow {X)) _ v we find 

(i* oup (A))_„ = »,'(r A - 1 (L"™'(A))J. 

The assertion follows from this and (6.30). □ 

Hereafter crystal lattice will mean the upper crystal lattice at q — (we drop the 
superscript). 

A basic property of the VOs (6.11a), (6.11c) is that they preserve the crystal lattice 
[4]. Fixing the lowest weight vector \u* x ) e V* a (X) such that (\u\},\u* x )) = 1 we normalize 
(6.11c) as 

$V>±®K» - l«A> + -"- 

Let $ £j „ be the weight components of $ = and similarly let Qy* (v± <E> •) = 

Proposition 6.3. Notations being as above, we have for all e and n 

(0 

$ e „L(A) C L{ji), (6.33a) 
$:„L* a ( M ) c L* a (X). (6.336) 

(ii) For some s > we have 

$ £ „(L(A) A _ e ) C q-^- 2 ^- s L(p) (6.34a) 
$*„(L* a ( M )_ M+l ) C j-s»-%f)-' £ '«(A) (6.346) 

Proof. Assertion (6.33a) is proved in [4]. The argument in [4] shows also that (6.34a) 
is true if it holds for £ = 0. In Appendix 3 we verify the latter statement explicitly. 

From the normalization of <&y^ we have the relation 

= (**)*> ** = (id ® C-)$(a 2 ) x (-a) or 1 
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according as A = A or Ai. We can verify further that 

(T- 1 <8> id) ofoT A = ±($+ <g> v*_ - $_ ® 

where we put $ = <I> + ® v + + $_ <g> w_. In view of (6.31) the assertions (6.336), (6.346) 
follow from these relations. □ 

Thanks to (ii), is well defined on the g-adic completion 

Um(L(/i)®L)/(g n L(/i)®L), 

i.e., when applied to an infinite sum u = ^2 n u n such that u n E q M "L(X), lim^oo M n 
— oo, then for any N there are only a finite number of non-zero terms modulo q N in 
Hence we can iterate the VOs of type (6.11a) finitely many times (we drop the 
symbols for completion) 

L(X) *^m®L 



— > L(X) ®L®L — ► . • ■ ■ 

As mentioned already (Sect. 2) we conjecture that after suitably normalizing VO the 
infinite iteration also makes sense. 

6.8 n-point functions To examine the behavior of the embedding (6.35) we need to 
study the n point correlators as n tends to oo. We give below what is known to us 
about the general n point functions. In this subsection we consider VOs of the type 

$(z) = ** v (s). 

Let n = 2m or 2m — 1, A = A and n, = A or Ai according to whether n is even or 
odd. We define the vector w n (z) G V® n by 

(UA|*(zi)-"$OZn)Ki) 

n m-l , g , s 

= n ^ iy/i n ^ ^ m,i n 00r * ^ , *o. ^ 

j = l j=l i<j ^ 31 

A priori w n (z) ]J 1 j l =i( Z2 3- lZ2 j)~ m+: ' * s ^ en a f° rma l power series in z n ,z n /z n -\, 
Z2/Z1. In general, let $ : V(p) — ► V(A) ® be an intcrtwiner (where W is a finite 

dimensional module) and set $ua = ma ® w + • • •. Then w must satisfy [4] ej A,/l ^ +1 w = 
for all i. In the present case W = V Z1 <g> • • • <g> V^ n , A = A . Because V has a perfect 
crystal, so does V®" [2] . This implies that up to scalar the linear equations 

A(e ) 2 w(z) - 0, A( ei )w(z) = 0. 

admit at most one solution (recall that the first equation involves the indeterminates Zi), 
and hence that our w{z) is a rational function in zi, ■ ■ ■ , z n up to an overall scalar factor. 

Set R(z) = PR(z) x (1 - q 2 z)/(z - q 2 ) with R(z) defined in (6.16). Then we have 
the following properties of w n {z): 



Ra+i(zi/z i+1 )w n (z) = (siw n )(z) i = 1,- • -,n - 1 



(6.37a) 
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where Sj = and we set (sf)(zi, • • • , z n ) = f{z s -i^, • • • , 2 s -i( n )) for a permutation 

s. 

w„(g 6 zi, z 2 , •••,z„) = Aig 3m - 3 Pi„ • • • Pi 3 P 12 w n (z 2 , ■ ■ ■ , z n , zx) (6.376) 
where A\ = -t^ 1 for n = 2m, = (qt^ 1 ) 3 ^ 2 for n = 2m - 1. 

W2m(zi, • ' ' ,^2m-l,0) =lU2m-l(zi, " ' ' , Z 2m -l) <S> «- ~ <?(/l .W2m-1 (^1 , ' ' ' j Z2ro-l)) ® V+, 

2m-2 

W2m-l(zi,- ■ ■ ,Z2 m -2,0) =W 2m -2(zi,- ■ ■ ,Z 2m -2) ®V + X ] [ Zj. 

(6.37c) 

Property (6.37a) is a consequence of the commutation relations of VO. Property (6.376) 
is a reduced form of the gKZ equation under (6.37a). Finally property (6.37c) follows 
from the normalization of the VOs (6.8a). 

Let us introduce the coefficients a n (s\z) by 

w n (zi, ■■■,z n ) = ^a n (e\z)v El <g> • •• ®v £n . 

e 

The e = (ei, • • • , £„) range over £j = ± such that the number of +'s equals m. Then 
(6.37a) is rewritten in the form 

o(- ■■-+■■■ \z) = C 7 i a(---+ |z) (6.38) 

o(- ••££•• • |z) = q<Tia(- ••££••• |z). 

Here the operator cr^ is given by 

One can verify that they obey the Hecke algebra relations (Lusztig [34]) 

a l o- i+1 o- i = o-j+iCTiCTj+i, o-jCTj = o-jO-j (\i - j\ > 1), (cri + q)(<Ji - q^ 1 ) = 0. 

Using (6.38) the coefficients a(£|z) arc determined uniquely once we know e.g. ao(z) = 
a(e\z) for £i = • • • = £ m = +, £ m +i = •••=£„ = —. From computations for small n we 
expect that w(z\, • • • , z n ) is actually a polynomial in and q, and that ao(z) is given by 

( _ 9) m(m-i)/2 _ ^ JJ _ ^.J x 1 fern = 2m, 

l<i<j<m m+l<i<j<n 

x z„ l+ i • • • z 2m _i for n = 2m - 1. 

(6.39) 

The representation of the Hecke algebra generated by this polynomial is irreducible for 
generic g; at g = 1 it specializes to the Specht module [35] of the symmetric group 
associated with the Young diagram of shape (m, m) or (m, m — 1). We hope to be able 
to extract exact information about the asymptotics of n point functions by analysing 
these representations. 
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§7. Space of states 

7.1 Vacuum, shift and energy In Sect. 1 through Sect. 5 we have emphasized the role 
of the space t^(A) <8> V(\')* a where A, A' G {Ao, Ai}. To treat the infinite sums of vectors 
safely, we introduce its g-adic completion. Define 

C x ,\i = limi(A) <g> L* a (\')/q n L(\) ® L* a (A'), 

Va : A' = Q({q)) ®Q[[g]] A,A'- 

Here Q((</)) (resp. Q [[<?]]) denotes the field (resp. ring) of formal Laurent (resp. power) 
series in q. The space Va,a' inherits the left {/-module structure of level 0. Let V A A , 

denote the right [/-module structure on Va,a' given by f.x = a~ 1 (x).f (/ G Va,a',^ G U). 
If we regard an element / G Va',a (resp. g G V A A /) as a linear map V(X) — ► V(A') (resp. 
V r (A) — > T^ r (A')), the left (resp. right) [/-action is given by the adjoint action: 

ad(x).f = ^X(i) o/oa(i( 2) ), 

g.ad r (x) = ^a _1 (a;(2)) ojoi (1) , 
where x G U. We define the bilinear pairing V A A / x Va',a — ► Q((<z)) by 

( 9 \f) = try(A)(g : 4pg 4 ; /) - (7.D 

It is non-degenerate and enjoys the property 

(g.ad r (x)\f) = (g\ad(x).f) Vx G U. (7.2) 
Note that the denominator in (7.1) 

g-(4p,A) 

1*V(A)(«- ') = M 

is the principally specialized character of V(A). 

For each v we take bases {\u\ v) }} C L(A)„, flu^*)} C (i(A)„)* which are dual 
with respect to the canonical pairing ( , ). We call the canonical element 

|vac) A = ^K M )®K M *) gVa.a (7.3) 

i,i> 

the vacuum. This corresponds to the identity map V"(A) — ► V{\). Hence it is clear that 
the vacuum gives rise to the trivial representation over U: 

ad(x).\v&c)\ = e(x)|vac)A Vx G U. 

Since (\u^ v) *}, T x \u { ? ] )) G q( 2 P> x ~ v )A we see from (6.31) that |vac) A in fact belongs to 
Va.a- Define the right vacuum A(vac| G V A A analogously. The pairing (7.1) is so normal- 
ized that A(vac|vac)A = 1. 

Consider the composition of the maps 

L(A) ® L* a (X')L(fi) ® L <g> L* a (X')L(fx) ® L* a (fi'). 

By Proposition 6.3 it extends to the map T : V\ a' — > V^'. Now let A' = A and fi' = fi. 
We define the translation and energy operators T, H by 

T = -gjT (7.4) 

H = -{q-q- x ){T odoT' 1 - d) (7.5) 

where d is given in (6.10). 

The following is the first property we expect for the vacuum vector: 
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Proposition 7.1. The vacuum vector is invariant under T, H: 

T\vac) x = |vac) M , H\vac) x = 0. 

Proof. The second equation is clear from the first and d\v&c)\ = 0. Put <& A y = E ® 
Ve , (*$J* = $^ y * a 1 = E ® «e- :t suffices to show that 

-gfi |vac) A = -^^$ e | U M) ® {** e M v >) 

coincides with Ivac)^ (note that (<&*)* sends V(A)* to V(fj,)*). This is equivalent to 
showing 

^ s \u^}m) t \u^*},\v)) = -(g,)-'\v) 
for any \v) € V(/z). Since the left hand side is 

=£Mk v) ><k ,0 *>.*>>> 

= £<M> e », 

the assertion follows from (6.26). □ 

7.2 Creation and annihilation operators Let us define the components of the VOs of 
type (6.12a), (6.126) by 

$ A M (z) = v+ ® $+(z) + ® $-(z), 
*±(*) = *£a(*)(»±®(0) 

Hence from (6.136) we have the relations 

$;(z) = -^.(V), (*) = $+{zq 2 ). 
Using these operators, we would like to define the n-particle states to be 

(zx)-.. $*>„)» id) |vac) A (7.6) 

with £i,---,£„ e {+, — } and \z\\ = ■ ■ ■ = \z n \ = 1. The n-particle states in the dual 
space are to be defined analogously, using \(vac\ and <b e (z) in place of |vac)A and $£(2) 
respectively. These are eigenstates of the shift and the energy operators (see 7.2). 

Let us examine the meaning of (7.6) more closely by studying the Fourier compo- 
nents 

dz\ dz n m 



| Zl |=...=|z n |=i 27Tizi 2iriz n 1 
x (<&*>)••• ^ (* n ) ® id) |vac) A , (7.7) 
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where the mi range over all integers. For dcfinitcncss we take n = 2. Eq. (7.7) is an 
infinite sum over i, v of the terms 

Now for each (u| e V r (A), g V(A) the function (zi)<I>* 2 (^2)!"') is convergent 

in the domain \z\ \ 3> \z 2 \, \q\ < 1, and has a meromorphic continuation with respect to 
2ij z 2 g (C) x with poles only at Z1/22 = q~ 2 ,q 2 ,q 6 , ■ ■ ■ (see Proposition 6.1). Hence the 
integration in (7.8) is meaningful. Notice that (7.8) is different from applying the weight 
components <& £lTO of $ £ (z) 

By the definition the latter is obtained by taking the contour in (7.8) to be \z\\ \z 2 \, 
and because of the pole at \z\/z 2 \ = q~ 2 the two expressions give different answers. 

The type II VOs do not preserve the crystal lattice, so the individual terms (7.8) 
comprise negative powers of q. However computations suggest (see Appendix 4) that 
when we sum over i, v the negative powers all disappear, getting thereby a well defined 
element of V\ \ . More generally it can be shown using the g-KZ equation that each matrix 
element {u^zi) ■ ■ ■ $* n (*„)K>, <u|$ ei (*i) • ■ ■ $ £n (^)K) (H e V r (X), \u') g V(A)) 
is holomorphic on \z\ \ = ■ ■ ■ = \z n \ = 1. Fixing n let T(n), ^(n) denote the span of the 
vectors (7.7) and its dual counterpart, respectively. We expect that 

(1) T{n) C V, r r (n) C V r . 

(2) With respect to the inner product (7.1) the spaces T r {m) and T{n) with m ^ n are 
orthogonal. 

(3) (7.1) restricted to T r (n) x T(n) is non-degenerate. 
Assuming these, we define the physical spaces of states by 

^ = Co^)> ^■ = e~ ^ r (n). 

They are analogous to the usual Fock space of free particles. 

Define the creation operator ip*(z) : T(n) — > T{n + 1) on JF, \z\ = 1, by 

^^)(^ 1 (^ 1 )...$* n (^)®id)|vac) A 

= ($ £ (^* i (^ 1 )---$* n (^)®id)|vac) A . (7.9) 

Both sides are to be understood in the sense of the Fourier coefficients as in (7.7). The 
annihilation operator (p(z) : J-{n) — > Tin — 1) is defined to be the adjoint of the dual 
counterpart T r {n) — ► T r {n + 1) of (7.9). Thus ip £ (z)\va,c) — by definition, and 

A (vac|<^ £l (zi) ■ • • <p £n {z n )ip* e , n (z' n ) ■ ■ ■ p*, (z 1 )|vac) A 

= tr v(A) (V 4 "<5 £1 ( Zl ) ■ ■ ■ $ £ „ (*„)$*, (4) . . . $*, (4)) /tlvw (q~^ . 

Because <p±(z) act as identity on the second component of £(A) <g> L* a (X') the com- 
mutation relations (6.20a~6.20c) as operators on V(X) can be readily translated. We 
have 

ip ei (z 1 )ip £2 (z 2 ) =Yip £ < 2 (z 2 )ip £ < i (z 1 )(R vv (z 1 /z2))l^, (7.10a) 
vtMXte) = E^(^)^ 1 («i)(flv.v.(«i/^))5J (7-106) 
ip £l {zi)tp* 2 {z 2 ) = Y<p*, 2 {z 2 )<p e [(zi)(Rvv*{zi/z2))lf £ l + g\S £l£2 S( Zl / z 2 ). (7.10c) 
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Here the delta function 8(z) = ^2 neZ z ™ ar i ses because of the pole of (6.20c) at z\ = z^. 
Thus the VOs provide us with a lattice realization of the Zamolodchikov algebra. 

Remark. In the limit q — ► 1, z = q~ 2)3 / l7r , the R matrix (6.18), (6.19) reduces to the S 
matrix of the su(2) invariant Thirring model [36] 

Rvv{(3 > ~ Wl-l)r(l) /? — 7ri ■ 

The commutation relations (7. 10a)- (7. 10c) become those for the Zamolodchikov opera- 
tors (up to rescaling the operators (fi e {z), ip^(z)). 

Having set up the mathematical definitions of the creation-annihilation operators 
we can discuss their transformation properties under the shift and the energy operators. 

Proposition 7.2. 

Ttp^T- 1 = T(z)<p±(z), T V * ± (z)T- 1 = -r(z)-ViW, (7.11) 

[H, <p±(z)} = -e(z)<p±(z), [H, <p* ± (z)} = e(z)<p* ± {z), (7.12) 

where 

«w-(«-«-')»5K«^)- 

Proof. We show (7.11), (7.12) for ip±(z). To see (7.11) we need to prove $^ A (z)T = 
r(z)~ 1 T^Y ft (z) , where the left and right hand sides are the compositions of 

L(A) ® L* a (A) — >L(ji) ®L® L* a {\)—>L{n) ® L* a (n)—>L <g> L(X) ® L* a ( M ) 

and 

L(X) ® L* a (A) — >L ® L(^) ® L* a (A) — >L ® L(A) ® L <8> L* a (A) — >L ® L(A) ® £* a (M) 

respectively. This follows from (6.24) by multiplying $y 2 x* (22) ® id from the left and 
setting z 2 = 1. 

Using (7.11) we calculate 

TodoT- 1 (p ± {z) = r(z)" 1 To (\d,ip±(z)] + <p±(z)d) or 1 

d_ 
dz 

The first term can be written as 

d_ 

dz 



-t(z)~ 1 T o (z— - A M + A x )<p±{z)) o T- 1 + tp±(z)T odo T- 1 . 



- t{z)-\z- - A M + A A )(T o ^ ± (z) o T- 1 ) 



= -(z-^\ogT{z))ip±{z) - (z-^- - A p + A x )<p±(z). 

The relations (7.12) follow from these. □ 

Let us compare them with the formulas for the energy and momentum of 'spin 
waves' obtained in refs. [10,37] via the Bethe Ansatz method. The result is given in 
terms of elliptic functions of nomc — q = e~ 7 as ([10], eqs.(26),(27)) 

e{6) = ^— sinh7dn ( -— #, k) , p(9) = am (^—9, k) — ^-ir. 



/ 1 7 1 7 r \ / 1 7 1 r, 

TT \ TT J \ TT J 2 

Identifying z — —e 2l ° and using the identity ([38], p. 509) 



e -i P (0) = 1/2 ^ 

Ml), (7.12) are precisely the sam 
favor of our mathematical framework of the particle picture. 



Q(qz 1 )' dz q — q 1 

we see that our (7.11), (7.12) are precisely the same. This gives a strong evidence in 
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§8. Discussions 

Let us summarize the content of this paper. 

1. We conjectured that the embedding [3] of the irreducible highest weight U q (st(2))- 
module V(Ao) into the semi-infinite tensor product, is given by the limit of n-point 
correlation functions of the U q (si(2)) vertex operators as n — > oo. We do not know the 
exact form of the correct normalization of the vertex operator. 

2. We postulated that the mathematical content of the infinite tensor product on which 
the XXZ Hamiltonian (after a suitable renormalization) is acting, is 

End c (^(Ao)e^(Ai)). 

For simplicity, we mainly used the half of this, i.e., Home (V(Ao), V(A ) © V(A\)). 

3. We made the dictionary between the physical and mathematical pictures, which 
contains: 

(a) The translation and energy operators (7.4), (7.5). 

(b) The inner product of states (7.1). 

(c) The vacuum states (7.3). 

(d) The creation and annihilation operators (7.9). 

4. By utilising the g-deformed Knizhnik-Zamolodchikov equation, we got the following. 

(e) The energy and momentum of the creation and annihilation operators (Proposition 
7.2). 

(f) The commutation relations of the creation and annihilation operators among them- 
selves, i.e., the Zamolodchikov algebra (7.10a)-(7.10c). 

(g) A conjectural formula for the n-point correlation functions of the q-deformed vertex 
operators for arbitrary n (6.36), (6.39). (This is not the same with the n-point 
correlation functions of the creation and annihilation operators.) 

5. We checked the validity of our picture on the following points. 

(h) A one-line proof of the character expression for the one-point function (in the sense 
of [23]) of the six-vertex model (Sect. 5). 

(i) Comparison of the energy and momentum of particles with the result obtained by 
the Bethe Ansatz (Sect. 6.7). 

(j) Comparison of the vacuum and two-particle states with the Bethe Ansatz eigenvec- 
tors in the anisotropic limit A = — oo (Appendix 5). 
(k) Comparison of the commutation relations of the creation and annihilation operators 
with the ^-matrix of the su(2)-Thirring model (Sect. 7.2). 

6. The following is a list of unsolved problems (besides the several conjectures stated 
in the main body of the paper). 

• For the XXZ model (and the six- vertex model), we want to know about the form 
factors of the local operators, the staggered polarization [27] and other correlation 
functions. 

Other models, for which a similar diagonalization scheme might apply, are, especially: 

• The vertex models with perfect crystals. The framework of the present paper admits 
immediate extension to the general case. Reshetikhin [39] proposes that the space of 
n-particle states in the generalized XXX model with higher spin is a tensor product 
of two factors, (C 2 )®" and the space of paths of length n of an RSOS model. It 
would be an interesting problem to examine this picture from our viewpoint. 

• The XYZ model, the hard hexagon model [13], Kashiwara-Miwa's model [40,41], 
the chiral Potts model [42,43]. These are massive models, each within a different 
category. No obvious extension of our scheme to any one of them is known. 
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Appendix 1. Action of U q (sl(2)) on V(A ) 

We will recall some basic properties of the global base of the irreducible highest weight 
£/ 9 (jj)-modulc V(X), which are given in [25,44]. Then, we will compute the actions of the 

Chevalley generators on some vectors of the upper global base in the case g = st(2) and 
A = A . 

Let -B(A) be the crystal for V(X). (See Section 2 of [2], especially Definition 2.2.3 
and Theorem 2.4.1.) We use the maps £j, ipi : B(X) — > N. (Sec Definition 2.2.3 and 
(2.2.15) in [2].) Set k(b) = e^b) + ipi{b). This is the length of the s[(2)-string of color 
i through b. For b G B(X) we denote by G up (b) G V(X) the upper global base vector 
corresponding to b, and by G low (b) G V^A) the lower global base vector corresponding 
to b. (In [26], G up is written as G (see Sect. 5) and G low is written as G\ (see Sect. 4).) 

We use the following properties. 

(1) Duality. 

Let ( , ) be the symmetric bilinear form on V(X) given in Sect. 6.7, (6.27), (6.28). 
We have 

(G up (b),G low (b')) =5 W . {Al.l) 

Therefore, if 

e t G up (b) = J2^(b',b)G up (b'), fiG up {b) = £ cj(6', b)G up {b'), 
v v 

then we have 

ei G low {b) = ^c' t {b,b r )G low {b r ), f l G low {b) = ^rc l {b : b')G low {b l ). 
v v 

(2) Crystalline action. 

We have 

e t G up {b) = [e^G^ihb) + Y, E wG(b'), 

v 

fiCr»(b) = [ Vi (b)}iG up (fib) + Y, F wG{b'), 

b> 

where E l b , b and F b , b are zero if h(b') > h(b). We also have 

lim — ^ = lim r -*£ r = 0. 

(3) Positivity. 

We have £$, 6 , Ffr b G ©£° =0 Z> [n], where [n] means (q n - q~ n )/(q - q' 1 ). 
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(1) is implicitly given in [25]. (2) follows from Proposition 5.3.1 and (5.3.8-10) in 
[25]. (3) is proved in [44] (Theorem 11.5). 

Now we restrict to the case g = st(2) and A = A . The crystal B(A ) is identified 
with the set of paths: p = {p(k)} k >i is called a path if and only if p(k) = (+) or (— ), 
and 

p{k) = (+) if k is even and k ^> 0, 
= (— ) if k is odd and k 3> 0. 

For a path p, we define its signature 




in such a way that 

h > ■ ■ ■ > l m > 0, 

p(k + 1) = p(k) if and only if k e {h, . . . , l m }- 
We call i7i the depth of the corresponding path. 

Example. If p = (■ • -p(4) p(3) p(2) p(l)) = (•••+ + + +), then I = ^2^. The 

depth of p is 3. 

In Tables 1 and 2 below, if I is the signature of p, we will write I to represent G up (p). 
It is quite convenient to associate a parity symbol 




with a signature / in such a way that 

Sfe = if Ik = fc mod 2, 
= iil k ^k mod 2. 



Namely, we represent ^ = 




In Table 1 we list the actions of and fi(i — 0,1) on the upper global base labeled 
by / such that m < 3. The symbol <j) means the highest weight vector. We abbreviate 
signatures by the differences to the signature of the operand. For example, 



/iD=[2](i)H+(5) 
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really means 



for all positive odd l\. 



We have derived these formulas inductively by using the properties (1—3). 

In the right-hand-sides of these formulas, we allow a parity symbol to have length 
greater than the actual depth of the accompanied signature. Namely, if a signature I 
is accompanied by a parity symbol of length m, then I — (lk)i<k<m can be such that 
l\ > ■ ■ ■ > lj > = lj + i = ■ ■ ■ = l m , and it is considered as {h)i<k<j- If a signature 
which breaks even this condition does appear, (e.g., l\ — I2 > 0), we understand that 
the corresponding G up (p) to be zero. 

There are some cases such that a parity symbol to have length smaller than the 
depth of the accompanied signature. In those cases we understand the corresponding 
terms are zero. 

In the left-hand-sides of the formulas, we always assume the length of a parity 
symbol is equal to the depth of the accompanied signature. 

Appendix 2. Vertex operator V z <8> V^(A ) — > V(A\) 

Let V(Ai)(i = 0, 1) be the irreducible highest weight J7 g (s[(2))-module with the highest 
weight Ao, and let V z be the 2 dimensional [/-module given in Sect. 6.2. 

Denote by V(Aj) the direct product of the weight spaces of V(Aj). There is a unique 
[/-linear 

$:V Z ®V(A Q ) — >V(Ai) 

such that ® ma ),«Ai) = lj where ( , ) is as in Appendix 1. (To be precise $ is 

a U ® Q(g)[2;,z _1 ]-linear map V z <g>V(A ) -> V'(Ai) ® Q(q)[z, z^ 1 } (see Sect. 6.4), but we 
shall not bother writing Q(q)[z, z -1 ].) We list 

p' 



for the first seven paths p whose signatures (see Appendix 1) are <f>, (1), (2), (3), 




For convenience, we define $ CT : V z ® V^(A ) — ► V'(Ao) by <I> cr = cr$, where a is 
the Q(g)-linear map V(Ao) — > ^(Ai) induced by the Dynkin diagram automorphism of 
[/q(s[(2)) which exchanges the colors and 1. The list in Table 2 is for The vertex 
operator $ : V z <g> V(A ) — ► ^(Ai) is given by 

$( u± ®v) = a($ CT ((±) (8>u)), 

and the vertex operator $ c : V z V(Ai) — > V^(A ) is given by 

$ c (w+ ® u) = z^ a (v^ ® c(w)), $ c (u_ <g> u) = $ CT (w+ ® cr(u)). (-42.1) 



/i(MD=[2](ii + i)H+ 



Of 






• 
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The formulas below are obtained inductively by using the result of Appendix 1. For 
example, 

$-(<;_ ®(1)Q =$°(v_ ®/ o 0) 



n=0 



= ]T z"(2n + 1){[2](1) H+ r — } - ^ V 1 £ z ^ 2n ) E 



n=0 



n=0 



n=l 



Appendix 3. Vertex Operator V(A ) — > V(Ai) ® V 

We follow the notation in Appendix 1 except that now we use <& for 



and set 



$ : V(Ao) -> V(Ai)®V 



$(u) = ® w+ + <&_(v) ® w_. 



We list $J = ct$± in Table 3, where a is the Dynkin diagram automorphism. We use 
the normalization ($_(«a ),mAi) = 1- Let us prove 



oo 

<I>l(0) = 5> 3 "(2n) B 



n=0 



(A3.1) 



from which the rest of formulas are inductively derived by using 

$£(/(>«) = (T^/i*^") + <*±,+*-(«)> 
#£(/it>) - q* 1 h® a ±{v) + 6 ± ,-*l(v). 

The proof goes as follows. We have 

(<S>l(e 1 v),v') ± = (^ ± (v)Jov')+5 ±i+ (^_(v),t v'). (43.2) 

By using these formulas for u = 0, we find that 

(*£(0),i/)=O if«'£lm/ 2 UlmA 2 . 

In general, the lower global base G low (6) (6 G B) has the property that £i{b) > k if and 
only if G low (&) £ Im/f. In our situation, from this follows that 



G low (p)) - if e (p) > 2 or £l (p) > 2 



(43.3) 
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Now, for a given p G Vo, take any sequence (i^, . . . ,ii) such that p — fi N ■ ■ ■ fi ± Po- 
Again, in general, the actions of and fi on the global base are given by 

ei G low (6) = [Ub) + l]iG^{eib) + Y,Fb' b G low (b>), 

v 

/ 4 G low (6) - [ £i (6) + l] 4 G low (/ 4 6) + ]T ^, b G low (6'), 



where Fy b and Ey b are the same as those in Appendix 1. In short, the actions create the 

correction terms (i.e., terms other than the combinatorial ones G low (eib) and G low (/i&)) 
of the string length greater than the combinatorial term. Therefore, we conclude that 

($^),G low (p)) =c($^),/ iw ---/ il 0) 

with some non zero c G Q(q), and further that 

($J(0),G low (p)) ^0 

if and only if 

P = /o(/i/o)Vfor + 
= (A/o)^ for - . 

Finally, by using (A3. 2) (with v = (j)) recursively, we get (A3.1). 

Let us compute io( n } (p l+n ,pi) — {G{4>)\ <&°_ o • • • o <&°_ \G{4>)) up to q 8 for an arbitrary 

n 

n by using Table 3. Let Vi (i = 1, . . . , 6) be the upper global base corresponding to the 
following symbols. 



vi v 2 v 3 v 4 v 5 v 6 

<t> (2)H (4)H (?) 

We can extract the following from Table 3. 

$1 (vi ) = «i + q 3 v 2 + q 6 v 3 mod <? 8 , 

<S>-(v 2 ) = (q-q 3 + q 5 )v 1 + {-q 2 + q A )v 2 + q 3 v 3 + q 2 v 4 mod q 5 , 
&-(v 3 ) = q 2 vi — qv^modq 2 , 

{Vi) = -q 3 vi + v 2 - qv 3 + q 2 v 6 mod q 3 , 
®-(vs) = v 3 modq\ 
®-{v6) = v 4 modq 1 . 

Here ll modq k,, is to be understood for those terms except for v\. The wi-terms are 
given by "modg' c+1 ". The reason we truncated the expansions at order less than q 8 
is as follows. For example, starting from v\ and applying $1 repeatedly one will get 

«2-terms only with power q 3 at the lowest. Then, unless one get a wi-term in the next 
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application of <f>^, we get at least one more power of q until one does finally reach a 
t>i-term. Therefore, modg 5 is enough (8 — 3 — f = 4) except for the Wi-term. 
Let us denote by (1, ji, . . . ,jk) a process 

$^ tj>^ $^ 
vi -» v h -> > -y Jfc . 

For example, if k — 1 we have two processes (12) and (13) that contribute to the final 
answer. Along with coefficients, the k = 1 process gives rise to q 3 (l2) + q e (13). We then 
proceed for a larger k step by step. Since we consider only up to q s , this terminates in 
finite steps. Picking up the A:-processes ending at jk = 1, we have 

fc = (1), 

k = 2 (q 4 -q 6 + q S )(121) + q 8 (l3l), 

k = 3 (-q 6 + 2g 8 )(1221) + g 8 (1231) - g 8 (1241), 

k = 4 g 8 (12221) + {q 6 - g 8 )(12421) - g 8 (12431) - g 8 (13421), 

k = 5 - g 8 ((122421) + (123421) + (124221)), 

k = 6 g 8 ((1242421) + (1243421) + 1246421)). 

Taking combinatorial factors into consideration, we get (for n > 5) 

wW(p„,Po) = 1+ (n - l)(g 4 - <? 6 + 2<? 8 ) + (n - 2)(-g 6 + 2 g 8 ) 

+ (n - 3)(g 6 - 2<? 8 ) + (n - 4)(-3<j 8 ) + (n - 5)(3g 8 ) 
(n-2)(n-3) 8 
2 9 

= 1 + (n - 1)<? 4 - nq 6 + "fr " 1 ) ^ modg 10 . 

For smaller values of n, we get 

^ (2) (po,Po) = 1 + 9 4 - q 6 + q S modq w , 
w (3) (pi,Po) = l + 2< ? 4 -3< ? 6 + 6< ? 8 mod<? 10 , 
w (4) (Po,Po) = l + 3 9 4 -4 g 6 + 9 9 8 mod 9 10 . 

The results agree with those in 6.8 obtained by solving the q-KZ equation. Similarly, we 
calculate 



L0 (n) ((- ■■ + - + - - +),p Q ) = -q + q 3 - nq 5 + 2nq 7 , 



(«) 



Thus we get 



((••• + - + + - -),Po) = -q + 2q 3 - (n + l)g 5 + (3n + 5)q 7 . 



<(■ ■■ + - + - - +),Po) = -q + q 3 - q 5 + q\ 

•• + - + + - -),Po) = -q + 2<f 3 - V + 7g 7 , 



which are in agreement with the result in [3] . 
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Appendix 4. Embedding of the vacuum, 1 and 2 particle states 

First we give a supplement to Sect. 2 about the embedding V(Ao)* a — ► V <S> V ® • • ■. 
Let us denote by u the automorphism of the algebra U q (sl(2)) given by 

u(ei) = fi, w(fi) = e u u(q h ) = q~ h . 

Let Vi and V2 be left {7 9 (sl(2)) -modules. A Q(<7)-linear isomorphism 77 : V\ — ► V2 is 
called w-compatible if and only if 

uj(x)i](v) — rj(xv) 
for a; G J7 g (st(2)) and u £ Vi. The map (6.32) in Sect. 6.7 

»y : V(Ao) - V(A a )* a 

is w-compatible. The map 

7} : V — ► V v± v T 

is also w-compatible. Suppose that 77^ : > V- (i = 1, 2, 3) are w-compatible, and that 
<& : Vi — ► V2 ® V3 is an intertwiner. Then the map 

& : V{ — >Vl®Vi m(v) ^ (V3 ®m) oPo^(v) {A4.1) 

is also an intertwiner. Here, P is the transposition of the tensor components. For a given 
path p £ Vi we set 

p* = (-p(l) -p(2) -p(3)---). 

The set 7?* = £ 7^} is naturally the crystal of V(K l )* a . Let $ : V(Ai) -> 

V{K.i-i) ®V be the vertex operator given in Sect. 4. Then 

$' : F(A,)* a -► V ® y(Ai_ 4 )* Q 
is also a vertex operator. From this we can define an embedding 

1' : V(A )* a -^V ®V ®V 
The embeddings 1 and 1' commute with the u-compatible maps r\\ 

V(Ao) ■■■(g)V(8)V 

F(A )* a V<S>V<S>--- 
From (6.32) and (.41. 1) we see that the vacuum state \vac) embedded in W is given 

by 

\ vac )=Yl J2 (-l) ht(A °-^k lp+Aol2 ^ p+ ^ 2 ^G up (p)) ®v°^G low (p)). 
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As we have noted in Sect. 2, we have 

*(G" P (P))U = |P>- 
For G low (p)(p) we must shift the q power; 

L(qW-M*G 1 ™(p))\ g=0 = \p). 

Note that 

^_ 1 ^t(Ao-M) (? |p+Ao| 2 -|p+M| 2 g -|Ao| 2 + |A'| 2 _ (_ ? )h*(A -M)_ 

Therefore, the term (-I^Ao-^Ip+AoIMp+^I^GW^)) ® w o t(G low (p)) contributes to 
the sum only with power q ht ( A o-v)_ L e t. us compute \vac) up to order g 3 . We need the 
following five vectors; 



V\ 


= G up (-- 


•+- + -) 


= G low (- • • + 


+ -), 


V2 


= G up (-- 


•+- + +) 


= G low (- • • + 


- + +), 


V3 


= G up (-- 


•+--+) 


= i J G-(... 


+ - - +) 


V 4 


= G up (-- 


•+ + -+) 


= i J G-(... 


+ + - +) 


v 5 


= G up (-- 


•+ - -) 


= G low (- • • + 





We have 

5 

\vac) = ^(-<7) mW t(«i) ®uj o b{vi) modg 4 , 
i=l 

where {m(i)}i<j<5 = {0, 1, 2, 3, 3}. 

In the following table we show the coefficients of paths in the expansion of b(vi) 
(1 < i < 5), up to the relevant order for each i. If we write • • • as a part of a path, 
we mean that the indicated part of the path is identical with the ground-state-path 
Pq. If we write (• • •), we mean that the indicated part may be void; otherwise it must 
not be void. We write ± or =p, if the indicated column of the path differs from the 
corresponding column of p . Therefore, for example, • • • ± =F • • • actually represents the 

paths (•••+- + + --),(•••+-- + + -),(•••++ + -),(•••- + + - + -), 

etc., simultaneously. 

Path Coefficient 

Vl 

1 

— h -q + q 3 

± T • • • -q + 2q 3 

± T ± T("0 2g 2 

±T"-±T(---) q 2 
±±TT(--0 -q 3 
± t ± t ± t(--0 Sq 3 
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v 2 



± T 


± T ■ ■ ■ 


1 -<-i \ 


o 3 


± T 


•■•± T 


••± T(-") 


-5 


+ 






l- g 2 


+ — 


+ 






± T 


• • + 




-q 


H 


H h 




5g 2 


± qz 


± T ■ ■ ■ 


+ 




± T 


••+ - 


+ 


2g 2 


± T 


...± T 




g 2 


+ + 






g 2 


- + 






1 










- + 


- + 




-3g 








-9 


H 


+ 




1 



''3 

V4 
V5 

1 

From this we get the expansion of \vac) up to order q 3 . The result reads as follows. 

Path Coefficient 

1 

2g 2 



± 


T 






-q + 


± 


T 


± =F • • • 




2q* 


± 


T 


■■±T- 




q 2 


± 


T 


± T ± 


T ■ ■ ■ 


-5q 3 


± 


± 


T T ■ ■ ■ 




-q 3 


± 


T 


■■± T 


± =F • • • 


-2q 3 


± 


T 


± =F • • • 


± =F ■ ■ ■ 


-2q 3 


± 


T 


■■± T- 


■■± T • • • 


-q 3 



This agrees with the perturbation expansion expansion (see [31]). 



The 1-particle state with the quasi momentum u (z = e m ) is an embedding of V z 

into V(Ai) ® V(A a )* a , or further into W. We denote the latter embedding by (in 

general, o z „,...,z! for the n-particle embedding). The quantum affine symmetry discussed 
in Sect. 1 assures that the vectors obtained by this embedding are doubly degenerate 
eigenvectors of "H and its higher order relatives. Because of the property of the dual 
modules discussed in Sect. 5 of [2], the existence of the embedding is equivalent to the 
existence of the vertex operator <£>(z) : V z ® T^(A ) — > V(Ai). 

Let us compute this embedding at q = (i.e., only the crystal) by using the data in 
Table 2. We discuss in details only on the w+-componcnt. First consider the top term 

in l { z\v+)] 

OC 

= $>-(2n)n®0 

= (•••- + - +1 + 

+ z 2 (--'+ - + 



- + ) + z( + + -| + - + ) 

-|+ - + -...) + ••• 
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So, we get X)^Lo z ™ [[2«]] from the top term. Because of the translational covariance, 
we expect the whole limit t[ 1 \v + ) to be J2 n ez z n [[2n]]. I n the lowest order in q, the 
term (-) ht ( A »-^q\p+ A «\ 2 -\p+^ 2 r]OL(G low (p)) gives rise to (-q)' lt ( A °-^(|p)). Therefore, 
among 6 more datum in Table 2, only p = ( h — — +) and (• • • — + — +) give non- 
zero contribution at q — 0, i.e., 2]] and z _2 [[— 4]], respectively. So, it is consistent. 
Similarly, the v_-component has the limit J2 neZ z n [[2n + 1]]. 

The two-particle state is V Z2 ® V Z1 embedded in V^(A ) ® V(A Q )* a or in W. This 

(2) 

embedding «4 2 , Zl is obtained by a successive application of the vertex operators. For 
example, the spin 1 component is 

(_)^(Ao-p) (Z Ip+AoI 2 -|p+mI 2 $z2 ( v+ g, $zi ( w+ G up (p))) ® 7? o t(G low (p)|A4.2) 

Let us use signatures and parity symbols to represent the paths and the upper global 
base vectors (see Appendix 1). If p = (p, we have 

® cj>) = <j> + zi(2) □+ z 2 (4) □+■■■€ V(Ai). 

Then, by using (A2.1) 

#'(22)(«+ ® <t>) = ^(1) □+ ^(3) □+■■■€ y(A ) modg, 
$'(z 2 )( W+ ®z 1 (2) Qee-^Q 



^1^2 



&(z 2 )(v + ® z 2 (4) □) ee -z 2 (3) □ - z 2 z 2 











(5) 




+^ 2 4 (I) 





e V(A ) modg, 



1 E V(A )modq. 



Therefore, the contribution to (A4.2) at q = from p = is 

z 2 [[l,0]] + z 2 [[3,0]] + --- 

- *i[[l, 0]] + zt4[[3, 2]] + Zl zi[[5, 2]] + ■ ■ ■ 

- z 2 [[3, 0]] - z 2 z 2 [[3, 2]] + z 2 z 2 3 [[5, 4]] + z 2 z 2 4 [[7, 4]] + ■ • 

So, by the translational covariance, we expect 

^(z^™ +1 -zr +1 4)[[2m+l,2n]] 



(.44.3) 



as the whole answer. The datum in Table 2 turn out to be consistent with this. In 
Appendix 5, the comparison of this result with the Bethe Ansatz calculation is given. 
Similarly, we obtain 

4 2 2 ) 2 >_ ® v.)\ q=0 = £ (zr'zT - zTz^Wm^n-l]}. (AAA) 

m>n 
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The case spin is equal to zero, is somewhat tricky because the limit q = can be 

(2) 

taken only after summing up certain series. Let us consider Li 2 , Zl {v~ <£>i>+). The first term 
is p = <j>, and we want to compute &'(z2)(v- <E> &(v+ <8> <fi)), or ^ a (z2)(v+ <£) z™(2n) 
(see (A2.1) and Table 2). A cumbersome (or rather interesting) feature here is that the 
terms proportional to <p contains negative powers in q. Keeping the whole terms that 
are proportional to <j>, and neglecting all other terms that vanishes at q = 0, we have 



<T(z 2 )( W+ ®0) = </> + ]Tz"(2n) 0, 



^{z 2 ){v+®z l {2) 0) 



zi 1 



z 2 q 2 1 - q 4 

OG 



*i(2)E 



<D-(z 2 )( U+ ®^ (4)H: 



n=2 

Z2<7 2 



1 - g 2 1 - q 6 



-«?(4)H-^ 2 



I 1 00 



n=3 



Extrapolating, we have 



1+ 



zi 1 



z 2 q 2 I — q' 



+ 



z\ 



1 - q 4 l-q i 



q 2 (z 2 - Zi) z 2 - ziq 4 z 2 - ziq 8 



z 2 <T 



Zl Z 2 



ziq 2 z 2 



ziq° 



Therefore, after summation the terms proportional to 4> vanish at q = 0. (In fact this is 
exact as shown in Sect. 6.5.) With this understood, we get 



4 2 2 ),>_ ® v + )\ q=Q = J2 - «)[[2m,2n]], 



(A4.5) 



m>n 



and similarly, 



9 =0 



E 



,m+l 



m n+1 



Zl Z 



)[[2m+ l,2n+ 1]]. (A4.6) 
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Appendix 5. Comparison with the Bethe Ansatz States 

In this appendix we calculate the Bethe vectors at q — and compare them with the 
results in Appendix 3. We assume the periodic boundary condition and hence that the 
length N of row to be even. Set N — 2n. Following [27] we label the standard basis 
vectors v £l • • • <S> v eN of V® N by the locations of v_: (x\, ■ ■ ■ ,xm), I < x\ < ■ ■ ■ < 
xm < N. We denote by \x) = \x\, ■ ■ ■ ,Xm) € V® N the corresponding vector with spin 
n — M. It is known that the vector 

v = ^2 f{xi, ■ ■ ■ ,x M )\xi, ■ ■ -,x M ) (A5.1) 

is an eigenvector of the XXZ Hamiltonian H if f(x\, ■ ■ ■ , xm) is of the following form: 

M 

f(x u ---,x M )= n E (^).^(fe))n(- F Ki))) x '" 1 ( m - 2 ) 

<?es M i<j<k<M i=i 

where 

sm(u — v) sm(u — p^y) 

A = , g=-e 7 , 7GR>o, 

and {uj}^L 1 is a solution of the following Bethe Ansatz equation (BAE) 

M 

FK) 2 " = - JJ S(u k , Uj ) for 1 < k < M, (A5.3) 

3=1 

where 

sin(w — v + i-f) 



S(u,v) = 



sin(w — v — ij) 



We calculated the form of Bethe vectors at q = in the cases of 2/i-particles with 
spin h and ft, — 1. Here the particle number is the same as the number of holes in [10] or 
[9]. As a consequence we find that the Bethe vectors at q = coincide with the q — > 
limit of eigenvectors calculated in App. 4. The correspondences are given by 

• the ground states (Example 1), 

• two particle states with spin 1 (Example 2) ((A4.3)), 

• two particle states with spin (Example 3) ((A4.6) and (A4.5)). 

(i) The case of 2/i-particles and spin h. 

Here we consider the case M = n — h for some h G Z>i. It can be shown that there 
is a solution {uj}^~^ = {Aj}""/ 1 of (A5.3) which has the expansion of the form 

oo 

e 2 ^ = Xf (1 + ^ \fq l ) for l<j<n-h (A5A) 
»=i 



Diagonalization of XXZ Hamiltonian 51 

where q = — e -7 . Substituting (A5A) into (A5.3) and comparing the coefficients of q° 
we have 

n — h 

( A (0)j-(n+fc) = (.jjn-h+l -Q A (0) for 1 < j < n _ /j. (A5.5) 

fe=l 

Proposition A5.1. Let {rj} 2 ^ U {Mj'}"=i be a partition of {0, 1, ■ • • , n + h — 1} such 
that n < • • • < r 2 h and iti < • • • < fJ, n -h- Let 9 be a real number which satisfies 

ir(— 2/iui + y^?_i r,) 7r 

= — ^ 7?" 1 mod-Z, (A5.6) 

n(n + ti) n 

-,<e<-n+ 2n ^X 1) . (A5.7) 
n + h 

Set 

X^=e i9 , Af=A< 0) a" l<j<n-h, 

where a = exp(^2i.). Then {A^}™"^ is a solution of (A5.5) and any solution of (A5.5) 
is obtained in this way. 

Remark For each {rj} 2 ^, if n is sufficiently large, there are 2(/ii + 1) 6*'s which satisfy 
(A5.6) and (A5.7). More exactly let #o be such a that satisfies (A5.Q) and — 7r < 9 < 
-TT+f L . Then o ,j = #o + v (° < i < 2 Mi + 1) satisfies (A5.7) if n is sufficiently large. 
By expanding F(Aj) and E(Aj, Afe) in g at g = we have 

F(A J ) = -(Af )- 1 (l + 0(g)), 
A (o) 

E(Aj ; Afc) = 0) fc g -i(i + 0( g )). 

Substituting these expressions to (A5.2) we have 
Proposition A5.2. The eigenvector v is given by 

v = const. ^(A^a"" 1 ) - ^ fc=1 kj ■ II Pfe * ' a ) 
{fe} i=i 

x 2(-l)^("+ fc - 1 )(Ai 0) )-( n - h ) I1 a- !Bl(2 ' l ' il+ ^=i r< Vi. ' ' ' >**-/>) modg, 



where const, is an overall factor, and 



P l (a)= [] («'-<*'')> 

-(n+/i-l)<i<j<0,ij'#; 

l/(fc,a)=det(a-("- r *) fe 0i<«<2/ l - 
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Here and x\ varies subject to the following conditions. 

> k± > ■ ■ ■ > k 2h >-{n + h- 1), 

1 < xi < b„- h + n + h+ 1 (< 2/1+1), 

{b 3 )U = {-J^i^M^VU 6i > ■ ■ ■ > bn-h- 

For these data {x 2 , ••• , x 2 h} 'are determined by 

xj = xi + 2(j — 1) + i for —ki — i + 2 < j < —k i+ i — i and < i < 2h , 

where we consider ko = 1 and k^h+i = n — h. 
Lemma. Pk+i{a) = —aPk(a). 
Set 

lj = 2(kj + n) — n + j — x\ for 1 < j < 2h. 

Now we shall take a limit n — > oo in such a way that a Tj (1 < j < 2h) are finite and 

(— A^°')" has a limit. By the Lemma, if k\ — k 2 h is sufficiently small compared with n, 
the ratio Pk i (a)/Pk 1 (a) goes to one in the limit n — > oo. Using this and setting 

Zj = a r > (1 < j < 2ft) 

we have in the n — > oo limit 

Proposition A5.3. Suppose \l\ — l 2 h\ -C oo and -C oo. Then in the limit q — > 

and n — > oo, the coefficients v(x) of v v(x)\x) is 

u(or) = ±det(^ J ')i< iJ <2/ l (J|2:j) 2 . 

i=i 

where x\ — 1 or 2 and the signature is same if x\ is same. 



Example 1 (h = 0, Ground states) In this case a = exp(^p) , u.j = j — 1 (1 < j < n), 
A^ 0) = -1 or -exp(^) and \f ] = X^oP' 1 (1 < j < n). 

(1) A<°> = -1 

u = — 11, 3, 5, • • -,2n- 1) + |2,4,6, • • ■ ,2n) mod g. 

(2) A< 0) =-exp(^) 

v = |1, 3, 5, • --,2n- 1) + |2,4,6, ••-,2n) mod g. 



Let us write the vector of the form xi, • • • ,im) by using the chain of + and — . 
Namely, for example, the vectors |1, 3, 5, • • • , 2n — 1) and |2, 4, 6, • • • , 2n) are 

1 2 3 4 • • • 2n - 3 2n - 2 2n - 1 2n 
- + - + ••• - + - + 

+ - + -••• + - + - . 
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Example 2 (h = 1, 2-particle states, spin 1) The coefficients of the limit vectors are 

v(l,x 2 ,---,x n - 1 ) = ±{z*z$' -zfzfr), 

v(2, x 2 , ■ ■ ■ ,x n -!) = ±(z^zt - z^zt)(ziz 2 )i. 



Here 



and 



> ki > k 2 > -n, 



h = 2fci + n + 1 — x\, 
l 2 = 2k 2 + n + 2- x x . 

The numbers l\ and l 2 have the following pictorial meaning. We shall number the place 
between i — 1 and i by n+ 1 — i. Then l\ and l 2 are the numbers between + and + which 
occurs twice in the above picture. 

- + - + •••-+ | +-•••-+ | +-•••- + 

h h 



(ii) The case of 2ft,-particles and spin h — 1 (h > 1). 

Suppose that the solution of BAE consists of real quasi-momenta and one 2-string. 
([10]). Let us denote them by {Aj}".T 1 _1 and {z,w} respectively. Then BAE takes the 
form 

n — h — 1 

F(X 3 ) 2n = -S(\ j ,z)S(\ j ,w) [] s ( X J^k), (45.8) 

fe=i 

n-h-l 

F(z) 2n = S(z,w) J] 5(«,A fc ), (45.9) 
fe=i 

71 — h — 1 

FW 2 " = S(»,z) H S(w,X k ). (45.10) 
fe=i 

It can be shown that there is a set of solutions of (45.8)-(45.10) which has the following 
expansions at q = 0. 

oo 

e 2iX i =\f\l + 'jr i \f ) q i ) for l<j<n-h-l, (45.11) 

i=l 

oo 

e 2 4Z =z (0) g ( 1+z (n+fc-l) 9 n+h-l + ^ Z^V), (45.12) 

oo 

e ^iw =z (0) q -l^ + w (n+h-l) q n+h-l + J2 w {k) q k ). (45.13) 

k— n+h 
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Substituting these expressions into (A5.8)-(A5.10) we have, at q = 

n — h — l 



(\f)- n - h+i = (-i) n -v°>) 2 n A i 0) ' ( A5 - w ) 

fe=i 

(z (0) ) -2(/ l +l) = ( "-Q 1 A (0) ) 2^ (A5J5) 

fe=l 

We can solve (A5.14) and (A5.15) easily. 

Proposition A5.4. Let {rj} 2 ti U {/^j™"-^ -1 be a partition of {0, 1, ■ • • , n + h- 2} suci 
that n < • ■ • < r 2 h and /zi < • ■ • < Un-fr-i- $r aiJ d #c he reai numbers which satisfy 

0r = — T7 n mod -^-Z, (A5.16) 

nh(n + h - 1) nh v ; 

- 7 r<^<- 7 r+ 2 ^ 1 + 1) . (A5.17) 
n + h — 1 

n(-2hri 1 + Y] 2 _, r,) 

c = -(n-l)0 r + — , r J mod ttZ, (A5.18) 

n + h — 1 

-7T < 6» c < 7T. 

Set 

A (°) =e *V 2 (0) =e i9 ej A (0) =A (0) a w-m l<j<n-h-l, 

where a = cxp(^ rT ). Then {A^}?"/ 1-1 U {z, w} is a solution of (A5.14)and (A5.15). 
Any set of solutions of (^45.14)and (A5.15) is obtained in this way. 
Lemma. Let f(xi, • • • , x n ^h+i) be defined by (A5.2). Then 

f(x U X n - h+1 ) = f ( Xl , x n _ h+1 )q^ n - h -^ n - h -^ + ( g i(n-fc-l)(n-fc-4)) . 

Here fo{xi, ■ ■ ■ , x n -h+i) = unless (x\, ■ ■ ■ , x n -h+i) is one of the following form. 

(1 ) There exists p > 1 which satisfies 

x p+ i = x p + 1, Xi+i — Xi > 2 ifi^p, 

where x n -h+i ^2nifp=l. 

(2) xi = 1, x 2 = 2, x n - h+ i = 2n, x i+1 -Xi>2ifi^l. 

(3) xi = 1, x n - h = 2n-l, x n - h +i = 2n, ar i+ i - Xi > 2 if i ^ n - h. 

Let us calculate fo(xi, ■ ■ ■ , x n _h+i) m the case of (1) in Lemma. Let us again use the 
± description and denote by be the ± chain corresponding to (x\, ■ ■ ■ , x n -h+i). 

Lemma. Suppose that (p{i)) 2 ™i is of the following form. 

p(2m,i — i — 2 + xi) = p(2m,i — i — 1 + X\) = + for 1 < i < I — 1, 
p(2mi - I - 1 + x\) = p{2mi - I + x\) = 

p{2m, L -i + x 1 )= p{2m l - i + 1 + x{) = + for I + 1 < % < 2h, 
2 < mi < ■ • ■ < rri2h < n + h — 1. 
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Set 

kj = —rrij + 1 for 1 < j < 2h and X\ = 1, 2. 
Tien p = mi — I and, up to constants independent of {xj}™^ 1 

a ^ 2^=i JJ Pfei (o)D(n, • • • , r 2fc |fci, • • • , k 2h ). 
i=i 

where 

D(r lt ■ ■ -,r 2 h\ki, ■ ■ -,k 2h ) = det(a ( - z -V ( - : >-V)i 7 t ku ... t k 2h ,jjir 1 ,-,r2 h - 
Let us define kj (1 < j < 2h) by 

kj — n + 2k j + j — x\ for 1 < j < I — 1, 

fc; = n + 2ki + I — 1 — x\ , 

%• = n + 2kj +j-2-xt for / + 1 < j < 2/i. 
By taking the limits n — » exo, g — > in a similar manner as in Proposition A5.3 we 

have 

Proposition A5.5. Suppose \k\ — k 2 h\ <^ 00 and |fci| oo. Then in the limit q — > 

and n — > oo, the coefficients v(x) of v = J2 v{x)\x) is 

2h 

k ,T-r , X i 1 

u(a:) = ±det(V)i<i,j<2h(]_]_2j) 2h ■ 

where X\ — 1 or 2 and the signature is same as far as X\ is same. 

Example 3 (h = 1, 2-particlc states, spin 0) The coefficients of the limit vectors are 

v{l,x 2 , ■ ■ -,x n -i) = ±(% fel 4 2 - 4 2 4") for I = 1, 
v(l,x 2 , ■ ■ ■ ,x n -i) = ±{z^z 2 2 - z'[ 2 z 2 1 )(ziz 2 )~^ for I = 2, 
v(2, x 2 , ■ ■ •,!„_!) = ±(^4 2 - z^z^)( Zl z 2 )i for Z = 1, 
v(2,x 2 ,---,x n - 1 ) = ±^z\ 2 -z\ 2 z k J) for I = 2. 



Here 
and 



> fci > fc 2 > -n, 



kj = 2k j + n — x\, for 1 = 1, 
= 2k j +n+l-x 1 , for I = 2. 
The ±-chain pictures are 

- + - + ••• + - J - + •••- + J + -•••- + 

fci k 2 

for I — 1, xi = 1 and 

- + - + •••- + J +-... + - J _ + ..._ + 

fci fc 2 

for I = 2, x\ = 1. 
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Table 1. 



= ei = O / O 0=(1)Q A0 = O 

eoD=(-l)H eoH=0 eiD = ei0 =(-!)□ 



/oD=0 /o0=(l)D /iD=[2](l)H+(i) T /i0= J 7 



eo 
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eo 
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eo 
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eo 
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ei 
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ei 




= (») 




ei 
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ei 



= [2] 



/o 



fo 



*T = PI 1 " 



A@-M(i)S +p| {(! 



2 

-1 



^ n = pi (; Pi 



eo 



[3] 



+ [2] -1 



-2 



"2 I □ 



Diagonalization of XXZ Hamiltonian 



■59 




GO 



B. Davies et al. 




Diagonalization of XXZ Hamiltonian 



61 



Table 2. 



n=0 

oo 

<S> a [y- <8> <p) = zn ( 2n + !) □ 

n=0 

°° /9r> 

*'K®(i)D = E* B (T 

$-( W _®(l)Q = -z- 1 g -V + E^(2n + 2)0+f;z"(' 2n i +1 ] P 

n=0 n=l ^ ' ' ' 



$-(, + ®(2)B^z- 1 ^-(2)H+^^ j | ( 2 n + 2 )H + ^z>f Z Y 1 )U+E 



[2] 



n=0 



[2]' 



71=1 



2n+ 1 



71=2 



2/i 



*-(«_ ® (2) H) = -*-%(l) □ + E ^T2T( 2 " + 3 ) □+ E o 

71=0 ^ J 71=1 ^ 

$>+ ® (3) Q = - 



/2n + l\r 



71 = 2 



*'(»-8(3)a = -^ 2 ^f*-^'ii](2)[I]-(l) @ + £^(2n + 4)H + g 2 » I |( 2 " 1 +3 



( J) -) = ^- 1 ?- 2 (i) □ - (3) □ + E 



n=l 



2n+ 1 



00 / 2n 

EW 2 

71=2 \ 1 



'2n + l' 

$<T ( w -® (t) ft = E x " I 2 



$> + ®(4)H) 



-1 9" 



[2P([3]-1)^ [2] 2 ([3]-l) 



Zi>(2)H-(4)H-ff K T 









• 


(?) 


• 


-(i) 
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f = n q2[3] (2r. + 4) m+ f z»* ( 2U + 



(2n + 






• 




• 







n=2 v 



n— 3 



2n 



00 ^"^2 



\- *V[3] ,„ 
^[2]»([3]-l) l2n 



n=0 



n=2 V 



*'K®(?) £ 



)=-* 



-2 9 



-6 r -l„-2[ol °° ra 2 

[2]([3]-l) (2n + 4)H 



[2]([3]-l) r [2]([3]-l) 



E^( 2n 2 +2 )^+E*"| ^ 

71=1 ^ ' 1 1 n=2 \ 1 



2n 



z~ 2 q- 6 



[2]([3]-l) v [2]([3]-l) 



^ [2]([3]-l) 



oo / 2n + 2 \ 



2n+ 1 

E*"l 3 

n=2 \ 1 
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*+ W = - E « 3n+1 ( 2n + !) □ $ - W = E « 3 "( 2 ") H 



n=0 



n=0 



oo oo ^ 1 

n (a) o = ^ - E « 3n+1 ( 2n + 2 ) H - E i 3n+2 i 

n=0 n=l ^ 

*-((!)□) =£^ i ( 2 r) ff 

n=l ^ ' 



— 1 00 3n 00 /o 1 \ I 

\ ((2) ED = ^j-(l) □ - E ^-(2" + 3) □ - E * 3n+1 ( 2 ) E 



n=l 



*-((2)H) = 



[2]" [2]' 



1 J V 2 )H+E^r(2- + 2)H 



^2n + l^ ^ 



+ E« 

n=2 



.'in 



2n 



+E^" 1 ( i 

*+((3) O = ^ + ^f(2) H - E If (2n + 4) 



-,3n 



9 3 /3 



[2] WW ^ [2] V 



~ ^3n+l / 2n + gX 

h [2] I i J 



n=2 



2n + 1 
3 



[2] 



q 3n ~ 2 (2n + 2\ 

^Tl i ) 



n=2 
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-(( i ) -) = (i)n-ff(3)n+Efl 8B 



71=1 



_ 1 (2n + 1 
2 



+ 5> 3 " +1 | 2 



n=2 



+ 



q 2 f 3 



[2]*([3]-l)' 



-E 



[2] W U ^ [2] 2 ([3]-l) 



[2] 2 ([3]-l) 
■(2n + 5)n 



^^n/2n + 3\ □ V,,,,.:/ h ' 
^! M V 2 Jd ^ 
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